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INTRODUCTION 


The  in-bore  motion  (balloting  or  transverse  motion)  of  a  projectile  affects 
the  accuracy,  performance,  and  structural  integrity  of  the  projectile,  the  gun 
tube,  and  associated  devices.  Consequently,  it  is  an  important  factor  which  the 
designer  of  a  projectile  or  gun  system  needs  to  consider. 

Extensive  investigations  of  the  in-bore  motion  have  been  conducted  (ref  1 
through  5),  Equations  of  motion  of  a  projectile  launched  in  a  gun  tube  have  been 
formulated  and  many  projectiles  have  been  analyzed.  However,  the  past  work  em¬ 
phasized  only  the  formulation  of  the  problems  and  the  computations  for  several 
projectiles,  while  little  attention  was  paid  to  the  documentation  of  the  computer 
programs  used.  To  facilitate  future  analysis  and  to  aid  the  user  of  the  computer 
program,  the  computer  programs  listed  below  will  be  documented: 

o  BALLOT .  This  program  analyzes  the  in-bore  motion  (balloting  or 
transverse  motion)  of  a  projectile  launched  in  a  stationary  gun  tube  with  con¬ 
stant  twist  rifling. 

o  BALL0T2 .  This  program  computes  the  in-bore  motion  of  a  projectile 
launched  in  a  stationary  smooth  gun  tube.  This  type  of  gun  tube  has  straight 
rifling  or  no  rifling. 

o  BALL0T3 .  This  program  analyzes  the  in-bore  motion  of  a  projectile 
launched  in  a  transversely  vibrating  smooth  gun  tube  with  the  breech  end  fixed. 
The  gun  tube  has  straight  rifling  or  no  rifling. 

o  BALLOT 4 .  This  program  calculates  the  in-bore  motion  of  a  projectile 
tube  with  the  breech  end  fixed.  The  gun  tube  has  constant  twist  rifling. 

This  report  deals  with  the  program  BALL0T2.  The  other  programs  are  to  be 
documented  separately  in  the  future  if  necessary.  BALL0T2  is  intended  to  facili¬ 
tate  the  investigation  of  the  in-bore  motion  of  a  projectile  launched  in  a  smooth 
gun  tube.  The  term  "smooth  gun  tube"  means  that  the  gun  tube  has  straight  rifling 
only  or  no  rifling  at  all. 

This  report  documents  the  main  equations,  the  program  symbols  or  codes,  the 
flowcharts,  the  preparation  of  the  input  data,  the  output  of  results,  and  the 
operation  of  the  program.  A  sample  computation  is  included  to  show  the  data 
preparation  and  the  results  in  tables  and  figures. 

This  program  Is  coded  in  the  FORTRAN  IV  computer  language.  It  will  run  on 
CDC  6000  series  computers  and  requires  approximately  120000g  of  core  storage. 

This  is  the  first  version  of  the  program.  Some  of  the  statements  are  not  in 
their  simplest  forms.  In  order  to  make  the  program  more  easily  understood, 
reader’s  comments  are  cordially  welcome  to  guide  the  future  refinement  of  this 
program. 
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THEORY  OF  ANALYSIS 


The  computer  program  BALL0T2  analyzes  the  in-bore  motion  of  a  projectile 
launched  in  a  smooth  stationary  gun  tube.  The  mathematical  modeling  of  the  pro¬ 
blem  is  based  on  the  plane  motion  approach.  The  details  of  the  mathematical 
modeling  are  reported  separately  (ref  6)  and  not  repeated  here  except  for  the 
following  descriptions  included  to  aid  the  user  to  understand  the  program. 

When  a  projectile  is  launched  in  a  smooth  stationary  gun  tube,  there  is  no 
axial  spin  of  the  projectile  or  the  spin  is  very  small  and  therefore  may  be  ig¬ 
nored.  The  forces  acting  on  the  projectile  are  in  the  plane  containing  the  gun 
tube  axis  and  the  c.g.  of  the  projectile,  except  the  gravity  force.  The  effect 
of  the  gravity  force  outside  this  plane  is  small  and  may  be  ignored.  Thus  the 
projectile  may  be  considered  to  move  in  this  plane,  and  the  plane  motion  approach 
is  permissible  to  model  the  motion  of  the  projectile. 

The  plane  of  the  motion  is  represented  by  the  x-y  plane  (fig.  1).  The  x- 
axis  is  directed  along  the  center  line  of  the  gun  tube  and  the  y-axis  is  perpen¬ 
dicular  to  the  gun  tube  axis  and  points  upward.  The  origin  of  this  coordinate 
system  is  located  at  the  initial  intersection  of  the  gun  tube  axis  and  the  raid- 
section  of  the  rotating  band  (or  obturator).  The  projectile’s  c.g.  is  in  the  x-y 
plane,  which  may  be  vertical  or  not,  depending  on  the  inclination,  a,  of  the  gun 
tube  and  the  angle,  B,  which  the  x-y  plane  makes  with  a  vertical  plane.  The  z- 
axis  is  determined  according  to  the  right-handed  Cartesian  coordinate  system 
convention. 


The  forces  acting  on  the  projectile  (fig.  2)  are  the  base  pressure,  the 
contact  forces  and  the  frictional  forces  due  to  the  projectile-gun  tube  interac¬ 
tion  at  both  the  rotating  band  and  the  bourrelet,  the  gravity  force,  and  the 
aerodynamic  forces  which  resist  the  motion  of  the  projectile.  The  gun  tube- 
rotating  band  contact  pressure  acting  on  the  rotating  band  is  considered  to  be  a 
uniform  contact  pressure,  P  ,  and  an  unbalanced  normal  force  resultant,  N  ,  act¬ 
ing  at  the  maximum  deformation  point  of  contact  of  the  rotating  band.  The  base 
pressure  is  assumed  to  act  uniformly  at  the  projectile  base  and  the  rear  side  of 
the  rotating  band.  The  gas  leakage  is  ignored.  There  is  a  resisting  moment  at 
the  rotating  band  which  resists  the  tranverse  rotation  of  the  projectile. 


The  equations  of  motion  for  the  translation  and  rotation  of  the  projectile 
are  derived  from  the  plane  dynamics  principles  and  the  theory  of  elasticity  as 
follows  (ref  6): 
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Axial  distance  of  projectile  c.g.  to  aerodynamic 
center  (a.c.)  where  aerodynamic  forces  act.  It  is 
not  equal  to  the  free-flight  value  and  it  should  be 
determined  experimentally  under  the  in-bore  motion 
conditions • 

Reference  length  for  aerodynamic  moment  computation 

Coefficient  of  aerodynamic  drag,  lift,  moment.  It 
is  not  equal  to  the  free-flight  value  and  it  should 
be  determined  experimentally  under  the  in-bore  mo¬ 
tion  conditions. 

Angular  spring  constant  for  transverse  rotation  of 
projectile  or  moment  to  turn  projectile  transversely 
about  rotating  band  per  unit  angle 

Gravitational  constant 

Axial  distance  of  projectile  c.g.  to  leading  edge  of 
bourrelet 

Transverse  moment  of  inertia  of  projectile  about 

c.g. 


Spring  constant  of  rotating  band,  bourrelet 

Axial  distance  of  projectile  c.g.  to  rotating  band 
mid-section 

Mass  of  projectile 

Normal  force  resultant  acting  at  rotating  band, 
bourrelet 

Base  firing  pressure  per  unit  area 

Gun  tube-rotating  band  uniform  contact  pressure  per 
unit  area 

Bourrelet  radius 

Bore  radius  of  gun  tube 

Axial  distance  of  projectile  c.g.  to  projectile  base 
plane 

Velocity  of  rotating  band,  bourrelet  contact  point 
at  maximum  deformation 


w 


Width  of  rotating  band 


X,  X,  X 
XT 

y 


e 


0, 


e 


Axial  displacement,  velocity,  acceleration  of  pro¬ 
jectile 

Travel  of  projectile  when  rotating  band  exits  muzzle 

Normal  displacement  of  projectile  c.g. 

Maximum  normal  deflection  at  rotating  band,  bourre- 
let,  positive  (negative)  for  lower  (upper)  side 
contact 

Projectile  c.g.  eccentricity  from  geometric  axis  of 
projectile;  positive  (negative)  when  c.g.  is  above 
(under)  geometrical  axis  as  shown  in  x-y  plane 

Yaw  angle,  acceleration  of  projectile 

Frictional  coefficient  at  rotating  band,  bourrelet 

Air  density  used  to  compute  air  resistance  or  aero¬ 
dynamic  force  inside  the  gun  tube 


Any  standard  numerical  integration  method  may  be  used  to  solve  these  equa¬ 
tions.  The  Euler’s  method  and  the  constant  acceleration  technique  are  used  in 
this  program.  The  basic  equations  used  are: 
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Displacement,  velocity,  acceleration 
Time  step  size 

Subscript  to  denote  time  step  number 


(2a) 

(2b) 


The  frictional  coefficients,  the  spring  constants  and  the  average  gun  tube- 
rotating  band  uniform  contact  pressure  may  be  obtained  from  the  experimental 
data.  However,  when  the  experimental  data  are  not  available,  approximate  equa¬ 
tions  may  be  used  to  compute  the  data.  The  equations  are  derived  from  the  prin¬ 
ciples  of  elasticity  and  thin  shell  theory  and  are  as  follows: 
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ri>  re 
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Young fs  modulus  of  bourrelet,  rotating  band  material 
Wall  thickness  of  bourrelet 

Inner,  outer  radius  of  rotating  band.  In  case  of  a 
rotating  band  which  is  fitted  on  the  projectile,  the 
inner  radius  is  equal  to  the  original  value  before 
fitting  on  the  projectile. 

Radius  of  rotating  band  seat  of  projectile 

Poisson’s  ratio  of  rotating  band  material 


The  computed  uniform  contact  pressure,  pQ,  is  used  for  the  starting  period 
of  projectile  motion.  Experiments  show  that  the  contact  pressure  decreases  as 
the  projectile  moves  toward  the  muzzle.  A  linear  relationship  between  this  pres¬ 
sure  and  the  distance  to  the  muzzle  is  therefore  assumed  in  this  program.  The 
ratio  of  this  pressure  at  the  muzzle  and  at  the  beginning  of  projectile  motion 
may  be  estimated  according  to  known  experimental  data  of  similar  gun  and  projec¬ 
tile  systems. 


The  base  firing  pressure  is  used  in  the  computation.  If  the  given  data  is 
the  chamber  pressure,  it  is  converted  to  base  pressure  according  to  the  following 
equation  (ref  7): 
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Chamber  pressure 
Base  pressure 
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C  =  Weight  of  charge 

Wp  =  Weight  of  projectile 

Y  =  Ratio  of  specific  heats  of  propellant  gases  (its 

typical  value  is  1.25  which  is  used  in  the  present 
program) 

To  estimate  an  intermediate  value  from  two  given  data  pairs,  a  linear  inter¬ 
polation  equation  is  used  in  this  program.  This  equation  is  as  follows: 


where 


y*  =  y. 


yn  +  1  '  yn 


*n 


+  1 


(*i  '  Xn) 


(6) 


V  *n  +  1 


yn» 

y± 


Vn  +  1 


Given  x  values  at  known  time  or  space  stations  n  and 
n  +  1,  respectively 

Given  y  values  at  known  time  or  space  stations  n  and 
n  +  1 ,  respectively 

Known  intermediate  value  of  y 

Intermediate  value  of  y  to  be  computed 


FEATURES  OF  PROGRAM 


The  BALL0T2  program  computes  the  in-bore  motion  of  a  projectile  and  the 
associated  forces  during  the  launch  period  when  the  projectile  is  inside  a  sta¬ 
tionary  smooth  gun  tube.  The  program  may  compute  the  effects  of  the  following 
parameters : 

1.  Eccentricity  of  projectile  c.g. 

2.  Axial  c.g.  location 

3.  Bourrelet  diameter 

4.  Bourrelet  wall  thickness 

5.  Bourrelet  band  width  (length  of  bourrelet  portion  having  diameter 
same  as  bourrelet  diameter) 

6.  Bourrelet  location 

7.  Rotating  band  diameter 

8.  Rotating  band  width 
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9.  Rotating  band  thickness 

10.  Rotating  band  location 

11.  Length  of  projectile 

12.  Weight  of  projectile 

13.  Initial  position  of  projectile 

14.  Bore  diameter 

15.  Travel  of  projectile  or  the  length  of  gun  tube 

16.  Inclination  of  gun  tube 

17.  Friction  at  the  bourrelet 

18.  Friction  at  the  rotating  band 

19.  Air  resistance 

20.  Magnitude  and  variation  of  firing  pressure 

21.  Material  properties  of  projectile 

22.  State  of  wear  of  gun  tube 

The  program  requires  the  input  data  related  to  these  parameters.  The  input 
data  may  be  in  either  English  or  metric  units  of  measurement.  However,  they  are 
converted  to  English  units  in  the  computation.  The  outputs  are  expressed  in  both 
English  and  metric  units. 

The  program  computes  the  rotating  band  gun  tube  contact  pressure  and  the 
spring  constants  at  the  rotating  band  and  the  bourrelet  when  no  experimental  data 
are  available. 

The  input  firing  pressure  may  be  either  base  or  chamber  pressure.  If  cham¬ 
ber  pressure  is  used,  it  is  converted  to  base  pressure  in  the  program.  The  input 
data  which  are  related  to  the  locations  of  the  projectile  in  the  gun  tube  are 
expressed  in  terms  of  distances  to  the  muzzle.  This  facilitates  the  preparation 
of  input  data  for  the  investigation  of  cases  including  similar  gun  tubes  of  dif¬ 
ferent  lengths.  These  distances  are  converted  to  those  referred  to  in  the  ini¬ 
tial  position  of  the  projectile  c.g.  in  the  computation. 

The  results  of  the  computations  are  arranged  to  be  printed  out  in  tabular 
form,  plotted  in  figures,  and  recorded  on  a  magnetic  tape. 


OUTLINE  OF  PROGRAM 


The  program  first  reads  and  then  prints  the  input  data  cards  in  original 
order  for  data  check  and  future  reference.  The  main  input  data  are  later  printed 
in  tabular  form  for  easy  reference.  During  the  computation,  some  of  the  main 
computed  parameters  are  printed  at  specific  time  intervals  to  provide  a  check  on 
the  state  of  the  computation. 

The  maximum  and  minimum  values  of  the  results  of  the  computation  are  printed 
in  tabular  form  and  the  detailed  outputs  of  all  the  main  computed  quantities  are 
recorded  on  magnetic  tape.  This  information  can  also  be  printed  in  tables  in 
order  of  increasing  time,  if  these  data  are  requested.  Some  of  the  computed 
results,  such  as  axial  and  normal  displacements,  velocities  and  accelerations, 
normal  or  lateral  forces  at  the  rotating  band  and  bourrelet,  normal  accelerations 
at  axial  point  F  (any  desired  axial  point  such  as  the  fuze  c.g.  location),  yaw 
angles,  and  air  resistance  (aerodynamic  drag  forces)  are  plotted  in  curves 
against  time  and  projectile  travel. 

The  program  consists  of  one  main  program  and  seventeen  subroutines.  The 
principal  functions  of  the  main  program  and  the  subroutines  are  as  follows: 

1.  Main  program  BALL0T2 — Control  of  input  data,  computation  of  projec¬ 
tile  motion  and  forces,  and  control  of  output  data  and  figures. 

2.  Subroutine  CL1PV — Determination  of  an  intermediate  value  of  a  vari¬ 
able  associated  with  a  linearly  increased  argument  according  to  the  linear  inter¬ 
polation  principle. 

3.  Subroutine  CMAXMIN — Determination  of  the  extreme  values  of  computed 
variables  and  associated  time  and  projectile  displacement. 

4.  Subroutine  CRIPV — Determination  of  an  intermediate  value  of  a  vari¬ 
able  associated  with  a  randomly  varied  argument  according  to  the  linear  interpo¬ 
lation  principle. 


5. 

6. 

format  and  in 
7 

and  magnitude 

8. 

per  line. 

9. 

10. 

11. 


Subroutine  DATAINP — Input  of  program  data. 

Subroutine  DIPTBL — Printout  of  input  data  in  original  punched  card 
tabular  form. 

Subroutine  DTARO — Arrangement  of  input  data  pairs  in  proper  order 
to  be  used  in  computation. 

Subroutine  DTATBL — Printout  of  data  in  tabular  form  with  two  data 

Subroutine  FIGCTRL — Control  of  output  figure  drawing. 

Subroutine  FGCV3 — Curve  plotting. 

Subroutine  FGDA3 — Preparation  of  data  for  curve  plotting. 
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12.  Subroutine  FGMN3 — Determination  of  maximum  and  minimum  values  of 
axes  scales  for  curve  plotting. 

13.  Subroutine  FGPRTCV — Printer  curve  plotting. 

14.  Subroutine  FGPRTSB — Determination  of  curve  symbols  for  curve  print¬ 
ing. 

15.  Subroutine  FGSCLIR — Computation  of  proper  ratio  for  figure  axes 

scales. 

16.  Subroutine  FGSTRT — Computation  of  figure  axis  scale  starting  values 
and  units  for  curve  plotting. 

17.  Subroutine  OPDTA — Printout  of  all  computed  results  in  order  of 
increasing  time. 

18.  Subroutine  OPTBL — Printout  of  maximum,  minimum,  and  muzzle  values 
of  computed  results. 

19.  Subroutines  SCALE,  AXIS,  LINE,  SYMBOL,  PLOT,  and  ASPECT  are  library 
plotting  subroutines.  They  are  contained  in  the  FORTRAN  extended  library  of  the 
computer.  Other  subroutines  not  mentioned  here  but  used  in  the  program  are  all 
library  subroutines. 

The  outline  flowchart  of  the  program  is  shown  in  figure  3. 

Detailed  descriptions  of  the  main  program  and  subroutines  are  presented  in 
the  following  sections.  Descriptions  and  definitions  of  program  symbols  or  codes 
of  the  main  variables  or  parameters  and  flowcharts  are  included  in  each  section. 
The  intermediate  or  dummy  symbols  used  to  make  the  program  easier  to  follow  and 
the  symbols  for  subscripts  or  control  transfer  are  not  discussed. 

The  transfer  of  data  among  the  main  program  and  subroutines  is  through 
labeled  common  blocks  which  are  listed  in  table  1. 

A  discussion  of  the  main  program  and  subroutines  follows.  Note  that  the 
various  units  shown  are  for  reference  only.  They  may  be  changed  depending  on  the 
input  data  specif ication  or  control.  (See  input  data  coding,  page  15.) 


Main  Program  BALL0T2 


Computation  is  performed  in  the  main  program  BALL0T2.  It  uses  different 
subroutines  to  obtain  input  data,  print  computed  results,  and  plot  output  fig¬ 
ures.  The  definitions  of  program  symbols  are  presented  in  table  2. 

The  flowchart  of  the  main  program  is  shown  in  figure  4. 
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Subroutine  CLIPV 


This  subroutine  is  used  to  find  the  intermediate  value,  YO,  of  the  variable, 
YD,  corresponding  to  the  known  value,  XO,  of  the  linearly  increased  argument  XD, 
based  on  the  principle  of  linear  interpolation.  There  are  ND  pairs  of  known  data 
of  XD  and  YD.  The  known  value,  XO,  of  XD  determines  the  control  number,  N,  which 
in  turn  determines  the  value  of  the  slope,  SLP,  of  the  XD- YD  curve  of  the  compu¬ 
tation.  The  value,  YO,  is  then  computed  using  the  values  of  XO,  XD,  and  YD.  The 
value  of  SLP  is  saved  or  recomputed  for  the  next  computation,  depending  on  the 
values  of  XO,  XD,  and  YD.  The  computation  is  performed  according  to  equation 
6.  The  program  symbols  and  their  descriptions  are  listed  in  table  3  and  the 
flowchart  is  shown  in  figure  5. 


Subroutine  CMAXMIN 


Subroutine  CMAXMIN  is  used  to  determine  the  maximum  and  minimum  values  of  a 
computed  variable,  X  for  example,  and  the  related  time  and  axial  projectile  dis¬ 
placement.  These  values  are  determined  by  comparing  all  computed  values  during 
the  computation  process. 

The  program  symbols  for  the  various  quantities  are  listed  in  table  4  and  the 
flowchart  is  shown  in  figure  6. 


Subroutine  CRIPV 


Subroutine  CRIPV  determines  the  intermediate  value,  YO,  of  YD  corresponding 
to  the  known  value,  XO,  of  the  randomly  varied  argument,  XD,  based  on  the  princi¬ 
ple  of  linear  interpolation.  It  is  similar  to  subroutine  CLIPV  except  that  the 
independent  variable  is  now  varied  randomly.  The  descriptions  and  definitions  of 
symbols  used  are  the  same  as  those  of  subroutine  CLIPV.  The  program  stops  the 
computation  when  the  known  intermediate  value,  XO,  is  larger  than  the  given  maxi¬ 
mum  value  of  XD.  An  adjustment  of  the  time  step  size  may  be  necessary  to  rerun 
the  program  again.  The  flowchart  of  subroutine  CRIPV  is  shown  in  figure  7. 


Subroutine  DATAINP 


Subroutine  DATAINP  reads,  writes,  and  converts  input  data.  Input  data  may 
be  in  either  English  or  metric  units;  however,  in  processing  the  computation, 
they  are  converted  and  processed  in  English  units.  The  output  data  are  in  both 
English  and  metric  units.  Details  of  input  data  preparation  are  described  in 
Input  Data  Coding  (page  13). 

Most  of  the  program  symbols  used  in  this  subroutine  are  the  same  as  those 
used  in  the  main  program  and  have  been  listed  in  table  2.  The  additional  program 
symbols  not  presented  in  table  2  are  in  table  5.  The  flowchart  of  this  subrou¬ 
tine  is  shown  in  figure  8. 
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Subroutine  DIPTBL 


This  subroutine  prints  the  input  and  initial  data  in  tabular  form.  First  it 
prints  the  title  of  the  table,  and  then  it  calls  subroutine  DTATBL  to  convert  the 
data  to  proper  units  and  print  two  data  per  line  of  the  table.  The  program  sym¬ 
bols  are  the  same  as  those  mentioned  in  the  main  program  and  subroutine  DATAINP 
except  for  a  few  additional  symbols.  These  symbols  are  easy  to  understand  and 
are  not  discussed  here.  The  flowchart  of  this  subroutine  is  shown  in  figure  9. 


Subroutine  DTARO 


Subroutine  DTARO  reverses  the  order  of  X,Y  data  pairs  and  changes  the  magni¬ 
tude  of  Y  and  the  reference  point  or  origin  of  X,  if  desired.  For  example,  the 
muzzle  distance-gun  tube  diameter  data  may  be  converted  to  the  breech  end  dis¬ 
tance-gun  tube  radius  data.  The  program  symbols  for  the  various  main  quantities 
are  listed  in  table  6;  the  flowchart  is  shown  in  figure  10. 


Subroutine  DTATBL 


Subroutine  DTATBL  prints  the  given  data  pairs  in  tabular  form  with  two  data 
per  line.  The  given  data  pairs  are  multiplied  with  proper  conversion  ratios  so 
they  may  be  printed  out  in  both  metric  and  English  units.  The  descriptions  and 
definitions  of  additional  program  symbols  are  presented  in  table  7.  The  flow¬ 
chart  is  shown  in  figure  11. 


Subroutine  FIGCTRL 


Subroutine  FIGCTRL  controls  the  plotting  of  the  computed  projectile  motion 
and  forces  against  time  and  projectile  travel,  respectively.  The  plotted  figures 
portray  the  axial  and  normal  displacements,  velocities,  and  accelerations;  the 
yaw  angles,  velocities,  and  accelerations;  the  normal  or  lateral  forces  at  the 
rotating  band  and  the  bourrelet,  respectively;  the  normal  accelerations  at  a 
special  axial  point,  say  point  F,  where  the  acceleration  information  is  required; 
and  the  aerodynamic  drag  forces.  All  figures  are  plotted  in  metric  units.  Con¬ 
trol  numbers  are  provided  to  plot  the  desired  figures  and  plot  the  figures 
against  time  or  projectile  travel,  respectively.  The  main  program  symbols  are 
the  same  as  those  mentioned  in  the  main  program  and  subroutine  DATAINP.  Addi¬ 
tional  program  symbols  are  listed  in  table  8  and  the  flowchart  is  shown  in  figure 
12. 


Subroutine  FGCV3 


Subroutine  FGCV3  plots  one  to  three  curves  in  one  graph  with  data  prepared 
from  subroutines  FGDA3  and  FGMN3.  It  also  writes  the  figure  titles,  subtitles, 
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and  labels.  Additional  program  symbols  are  listed  in  table  9.  The  flowchart  is 
shown  in  figure  13. 


Subroutine  FGDA3 


This  subroutine  reads  and  converts  the  tape  recorded  data  to  suitable  values 
for  the  plotting  of  one  to  three  curves.  The  data  tape  is  labeled  tape  1.  Addi¬ 
tional  program  symbols  are  listed  in  table  10,  and  the  flowchart  is  shown  in 
figure  14. 


Subroutine  FGMN3 


This  subroutine  determines  the  maximum  and  minimum  values  of  axis  scales  for 
the  plotting  of  one  to  three  curves  in  one  graph  and  prepares  related  axis  la¬ 
bels.  These  values  are  determined  according  to  input  data  or  the  computed  re¬ 
sults  in  the  program.  Additional  program  symbols  used  in  this  subroutine  are 
listed  in  table  11,  and  the  flowchart  is  shown  in  figure  15. 


Subroutine  FGPRTCV 


Subroutine  FGPRTCV  prints  one  to  three  curves  on  output  paper  and  labels 
them  with  symbols  A,  B,  and  C,  respectively.  The  figure  fits  the  8  1/2  x  11  in. 
frame  size  and  is  suitable  to  be  included  in  a  report.  The  additional  program 
symbols  are  listed  in  table  12.  The  simplified  flowchart  is  shown  in  figure  16. 


Subroutine  FGPRTSB 


Subroutine  FGPRTSB  determines  curve  symbols  for  printer  plotting.  The  sym¬ 
bols  for  one  to  three  curves  are  prepared.  Additional  program  symbols  are  listed 
in  table  13  and  the  simplified  flowchart  is  shown  in  figure  17. 


Subroutine  FGSCLIR 


Subroutine  FGSCLIR  determines  a  proper  integer  ratio  for  figure  scales  when 
two  or  three  curves  are  plotted  in  one  graph.  However,  the  initial  computation 
of  this  ratio  is  first  done  in  subroutine  FGMN3  and  subroutine  FGSCLIR  refines 
this  ratio  when  it  is  larger  than  5.  The  symbol  N  is  first  used  to  represent  the 
originally  computed  scale  ratio  in  interger  form,  and  later  it  is  assigned  to  the 
recomputed  ratio.  The  symbol  ND  is  the  originally  computed  scale  in  floating 
point  form.  The  flowchart  of  this  subroutine  is  shown  in  figure  18. 
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Subroutine  FGSTRT 


Subroutine  FGSTRT  computes  the  starting  value  and  units  of  figure  axis  for 
curve  plotting.  It  determines  these  values  according  to  given  maximum  and  mini¬ 
mum  scale  values  or  maximum  and  minimum  values  of  computed  results.  The  program 
symbols  used  in  this  subroutine  are  listed  in  table  14  and  the  flowchart  is  shown 
in  figure  19. 


Subroutine  OPDTA 


This  subroutine  prints  all  computed  data  in  metric  units.  It  prints  these 
data  in  order  of  time  and  in  three  tables.  The  maximum  number  of  recorded  data 
is  600.  Less  than  600  data  points  may  be  printed  when  the  value  of  NPOPDTA  is 
greater  than  one.  The  printer  will  print  every  (NPOPDTA) th  data  from  about  600 
data  points  recorded  on  tape  1.  Additional  program  symbols  used  in  the  subrou¬ 
tine  are  listed  in  table  16;  the  flowchart  is  shown  in  figure  20. 


Subroutine  OPTBL 


Subroutine  OPTBL  prints  a  table  of  maximum,  minimum,  and  muzzle  values  of 
computed  results.  The  data  are  expressed  in  both  metric  and  English  units. 
Additional  program  symbols  are  used  to  represent  the  computed  results  in  differ¬ 
ent  units,  but  because  they  are  easy  to  understand  in  this  program,  they  are  not 
discussed  here.  The  flow  chart  is  shown  in  figure  21. 


PROGRAM  OUTPUT 


The  functions  of  the  main  program  and  various  subroutines  in  printing  the 
output  and  plotting  the  curves  of  computed  results  have  been  mentioned  in  the 
previous  sections.  The  following  presentation  is  provided  to  give  the  user  or 
reader  a  complete  idea  of  the  arrangement  of  the  program  output. 

The  program  first  prints  out  the  cover  page  of  the  output  sheets  and  the 
title  and  subtitle  of  the  analysis,  and  then  the  date  of  the  computer  run  is 
recorded. 

The  input  data  deck  is  printed  out  card  by  card  in  original  order  so  that  it 
is  easy  to  check  the  input  data  if  such  a  check  is  necessary. 

The  initial  projectile  position,  given  initial  yaw  angle,  and  computed  yaw 
angle  when  the  bourrelet  just  contacts  the  gun  tube,  are  recorded.  The  computed 
gun  tube-rotating  band  contact  pressure  will  be  printed  out  if  the  friction  coef¬ 
ficients  are  given  and  there  is  no  pressure  input  data. 
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The  dimensions  and  physical  properties  of  the  projectile  and  gun  tube  are- 
printed  out  in  a  table,  and  so  is  the  firing  pressure-time  data.  Similar  tables 
for  spring  constants  or  compression  test  data,  contact  pressure,  tube  wear,  and 
aerodynamic  coefficients  or  test  data  are  also  printed  if  there  are  input  data. 

When  the  program  starts  the  computation,  the  heading:  "Intermediate  print¬ 
out  of  computed  values  of  main  quantities  vs  time"  is  printed,  and  the  so-called 
main  quantities  are  printed  at  specified  time  intervals  (determined  by  input 
data).  It  is  useful  to  check  these  results  when  the  computation  fails  to  see 
which  parameters  are  abnormal.  The  total  number  of  computation  cycles  and  data 
recorded  on  magnetic  tape  are  recorded  at  the  end  of  this  table. 

At  the  end  of  the  computation,  maximum,  minimum,  and  muzzle  values  of  the 
main  computed  results  such  as  linear  and  angular  displacement,  velocity  and  ac¬ 
celeration,  and  forces  are  printed  out  in  a  table. 

Approximately  600  data  points  of  each  computed  main  quantity  are  recorded  on 
a  magnetic  tape  and  may  be  printed  out  in  three  tables  if  desired. 

Finally,  the  axial  and  normal  displacement,  velocity  and  acceleration,  yaw 
angle  and  velocity,  yaw  acceleration,  normal  acceleration  at  axial  point  F,  lat¬ 
eral  forces  at  the  rotating  band  and  the  bourrelet,  and  aerodynamic  drag  force 
curves  are  printed  out  against  time,  or  travel,  or  both.  They  are  also  plotted 
on  plain  or  grided  graph  paper  on  a  Calcomp  plotter.  The  printer  plotting  is 
printed  out  with  one  figure  per  page  of  computer  output  paper. 


INPUT  DATA  CODING 


Subroutine  DTAINP  reads  and  prints  input  data  in  original  order.  The  input 
data  may  be  expressed  in  either  English  or  metric  units.  However,  in  the  solu¬ 
tion  of  the  problem,  they  are  converted  and  processed  in  English  units  and  the 
output  data  are  in  both  English  and  metric  units.  The  following  input  data  cards 
are  necessary  to  define  a  problem  to  be  solved  by  the  RALL0T2  program.  The  sim¬ 
plified  diagram  of  a  projectile  (fig.  22)  indicates  the  needed  dimensions  except 
those  of  the  rotating  band  which  are  needed  for  average  gun  tube-rotating  band 
contact  pressure  computations.  The  input  data  cards  are  arranged  in  the  follow¬ 
ing  order: 


1.  Title  card  (1  card) 


List  TITLEM 

Format  (5A10) 
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Columns 

Symbol 

Format 

Description 

1-50 

TITLEM 

5A10 

Title  for  analysis 

2.  Subtitle  card  (1  card) 

List 

TITLES 

Format 

( 5A10) 

Columns 

Symbol 

Format 

Description 

1-50 

TITLES 

5A10 

Supplementary  information  about 
this  analysis;  if  no  information, 
leave  it  blank. 

3,  Computation 

and  output 

control 

card  (1  card) 

tlst 

NIDTBL ,  NCCK,  NOPTBL ,  NOOPDTA,  NP0PDTA,  NOAMS,  NDARAD 

Format 

(7110) 

Columns 

Symbol 

Format 

Description 

1— 10 

NIDTBL 

110 

Printout  of  input  data  in  tabular 
form;  if  =  0  or  blank,  printout 
input  data;  if  =  1 ,  no  printout. 

1 1-20 

NCCK 

110 

Printout  of  some  computed  values 
during  processing;  every  (NCCK)th 
data  recorded  on  magnetic  tape  is 
printed  out  for  intermediate 

-check  of  computations;  if  =  0  or 
blank,  NCCK  is  set  to  be  15  in 
the  program. 

21-50 

NOPTBL 

110 

Printout  of  maximum,  minimum,  and 
muzzle  values  of  computed  results 
in  tabular  form;  if  =  0  or  blank, 
printout  table;  if  “  1 ,  no  print¬ 
out  . 

5 1  -40 

NOOPDTA 

no 

Number  of  computed  data  to  be 

recorded  on  tape;  every 
(NOOPDTA) th  data  computed  will  be 
recorded.  If  =.  0  or  blank, 

NOOPDTA  is  set  to  be  such  a  num¬ 
ber  that  approximately  600  data 
are  recorded. 
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41-50 


NPOPDTA 


110 


51-60  NDAMS  110 


61-70  NDARAD  110 


Number  of  detailed  output  to  be 
printed;  every  ( NPOPDTA)  th  data 
recorded  on  tape  1  will  be  print¬ 
ed  out.  If  =  0  or  blank,  no 

printout . 

Units  used  in  projectile  and  gun 
tube  data;  if  =  0  or  blank,  data 
are  expressed  in  English  units; 
if  =  1,  metric  units. 

Unit  for  angles;  if  =  0  or  blank, 
angle  is  expressed  in  degrees;  if 
=  1 ,  radians . 


4,  Time  card  (1  card) 


List  TMAX,  TSEP,  TSLAST 

Format  (3F10.5) 


Columns 

1-10 


11-20 

21-30 


Symbol 

TMAX 


TSTEP 

TSLAST 


Format 

F10.5 


F10.5 

F10.5 


Description 


Estimated  time  duration  of  compu¬ 
tation,  ms;  computation  starts  at 
TIME  -  0  and  stops  at  TIME  =  TMAX 
or  when  projectile  is  out  of 
muzzle;  default  value  =  maximum 
time  of  firing  pressure-time 
input  data. 

Time  step  size  for  constant  step 
size  computation  or  first  step 
size  for  uniformly  decreased  step 
size  computation,  ms. 

Last  step  size  of  uniformly  de¬ 
creased  step  size  computation, 
ms . 


5.  Projectile  card  set  (3  cards) 
Card  A  Dimension  card 


List 


PR0JL,  RBX,  RBW,  RBD,  CGX,  CGD ,  CGEL,  CGEMT 


Format  (8F10.5) 


Columns  Symbol  Format 


Description 


PR0JL  F10.5  Length  of  projectile;  in  inches 

if  NDAMS  =  0  or  blank;  cm  if 

NDAMS  =  1. 
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11-20 

21-30 

31-40 


41-50 

51-60 


61-70 


71-80 


Card  B 

List 

Format 


RBX 


RBW 


RBD 


CGX 


CGD 


F10.5  Distance  from  projectile  base  to 

nearest  side  of  rotating  band;  in 
inches  if  NDAMS  =  0  or  blank;  cm 
if  NDAMS  =1. 

F10.5  Width  of  rotating  band;  in  inches 

if  NDAMS  =  0  or  blank;  cm  if 

NDAMS  =1. 

F10.5  Outside  diameter  of  compressed 

rotating  band  after  projectile  is 
loaded  in  gun  tube;  it  is  equal 
to  or  less  than  bore  diameter;  in 
inches  if  NDAMS  =  0  or  blank;  cm 
if  NDAMS  =1;  if  =  0  or  blank,  it 
is  set  to  be  equal  to  bore  diame¬ 
ter  in  the  program. 

F10.5  Distance  from  projectile  base  to 

c.g. ;  in  inches  if  NDAMS  =  0  or 
blank;  cm  if  NDAMS  =  1. 

F10.5  Projectile  diameter  at  c.g.  sec¬ 

tion;  in  inches  if  NDAMS  =  0  or 
blank;  cm  if  NDAMS  =1.  If  left 
blank,  it  is  set  to  be  equal  to 
bourrelet  diameter  in  the  pro¬ 
gram. 


CGEL  F10.5  Projectile  c.g.  eccentricity  in 

distance  from  c.g.  to  projectile 
axis;  in  inches  if  NDAMS  =  0  or 
blank;  cm  if  NDAMS  =  1 ;  if  left 

blank,  it  will  be  computed  from 
CGEMT  in  the  program. 

CGEMT  F10.5  Projectile  c.g.  eccentricity  in 

unbalanced  moment  about  projec¬ 
tile  axis;  in  in-oz  if  NDAMS  =  0 
or  blank;  cm-N  if  NDAMS  =  1;  if 
left  blank,  it  will  be  computed 
from  CGEL  in  the  program. 

Dimension  card  (continued) 

BRTX,  BRTW,  BRT0D,  BRTID,  BRTTK,  APAX,  APFX 
(7F10.5) 
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Columns 

Symbol 

Format 

1-10 

BRTX 

F10.5 

11-20 

BRTW 

F10.5 

21-30 

BRTOD 

F10.5 

31-40 

BRTID 

F10.5 

41-50 

BRTTK 

F10.5 

51-60 

APAX 

F10.5 

61-70 

APFX 

F10.5 

Card  C 

Physical  property  card 

List 

PWT,  TMOI 

Format 

(2F10.4) 

Columns 

Symbol 

Format 

1-10 

PWT 

F10.4 

Description 

Distance  from  projectile  base  to 
bourrelet;  in  inches  if  NDAMS  =  0 
or  blank;  cm  if  NDAMS  =  1. 

Width  of  bourrelet  portion  having 
same  diameter  as  bourrelet;  in 
inches  if  NDAMS  =  0  or  blank;  cm 
if  NDAMS  =  0  or  blank. 

Outside  diameter  of  bourrelet;  in 
inches  if  NDAMS  =  0  or  blank;’  cm 
if  NDAMS  =  1. 

Inside  diameter  of  bourrelet;  in 
inches  if  NDAMS  =  0  or  blank’;  cm 
if  NDAMS  =1.  if  left  blank,  it 
will  be  computed  from  outside 
bourrelet  diameter  and  bourrelet 
thickness  data. 

Thickness  of  projectile  wall  at 
bourrelet;  in  inches  if  NDAMS  =  0 
or  blank;  cm  if  NDAMS  =1.  if 
left  blank,  it  will  be  computed 
from  outside  and  inside  bourrelet 
diameter  values  in  the  program. 

Distance  from  projectile  nose  to 
projectile  axial  point  A,  where 
normal  acceleration  computation 
is  desired;  in  inches  if  NDAMS  = 
0  or  blank;  cm  if  NDAMS  =  1. 

Distance  from  projectile  nose  to 
projectile  axial  point  F,  where 
normal  acceleration  computation 
is  desired;  in  inches  if  NDAMS  = 

0  or  blank;  cm  if  NDAMS  =  1. 


Description 

Weight  of  projectile  in  lb  if 
NDAMS  =  0  or  blank;  kg  if  NDAMS  = 

1. 
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11-20 


TM0I 


F10.4 


Transverse  moment  of  inertia  of 
projectile  about  center  of  grav¬ 
ity;  in  lb-in^  if  NDAMS  =»  0  or 
blank;  kg-m2  if  NDAMS  *  1. 


6.  Gun  tube  card  (1  card) 

List 

ALPHA,  BORED,  FPEA, 

TRAVEL 

Format 

(4F10.5) 

Columns 

Symbol 

Format 

Description 

1-10 

ALPHA 

F10.3 

Gun  tube  inclination;  in  deg  if 
NDRAD  =  0  or  blank;  rad  if  NDRAD 
=  1. 

11-20 

BORED 

F10.5 

Bore  diameter;  in  inches  if  NDAMS 
=  0  or  blank;  cm  if  NDAMS  =  1. 

21-30 

FPEA 

F10.5 

Effective  acting  area  of  base 
pressure;  it  is  used  to  multiply 
base  pressure  to  get  pressure 
resultant*  in  in^  if  NDAMS  =  0  or 
blank;  cm^  if  NDAMS  *  1 . 

31-40 

TRAVEL 

F10.5 

Travel  of  projectile  from  its 
initial  position  to  rotating  band 
exit;  in  inches  if  NDAMS  =  0  or 
blank;  cm  if  NDAMS  «  1. 

7.  Projectile 

inclination 

card  ( 1  c 

ard) 

List 

BETA,  THETA,  THETAR 

Format 

(3F10.5) 

Columns 

Syrabo 1 

Format 

Description 

1-10 

BETA 

F10.5 

Angle  between  vertical  plane  and 
X-Y  plane  (or  plane  containing 
projectile  c.g.  and  gun  tube 
axis);  in  deg  if  NDRAD  =  0  or 

blank;  RAD  if  NDRAD  =  1. 

11-20 

THETA 

F10.5 

Initial  yaw  angle  (or  angle  be- 

tween  projectile  axis  and  gun 
tube  axis);  in  deg  if  NDRAD  =  0 
or  blank;  rad  if  NDRAD  =1. 
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THETAR  F10.5  Computation  of  yaw  angle  with 

bourrelet  just  contacting  gun 
tube  at  6  o'clock  position;  if  = 
1,  computation  will  be  done  in 
the  program  and  computed  value 
will  be  used  as  initial  yaw  angle 
in  the  computation;  if  =  0  or 

blank,  computation  is  also  done 
but  commputed  value  will  not  be 
used  as  initial  yaw  angle. 


8.  Figure  card  set  (1  to  13  cards) 

Card  A  Figure  plotting  control  card 


List 

NCGDVAA, 

NFGVSX 

NCGDVAN , 

Format 

(8110) 

Columns 

Symbol 

Format 

1-10 

NCGDVAA 

110 

11-20 

NCGDVAN 

no 

21-30 

NLF 

no 

31-40 

NAPFA 

no 

41-50 

NAF 

no 

51-60 

NYAWDV 

no 

NLF,  NAPFA,  NAF,  NYAWDV,  NYAWA , 


Description 

Plot  curves  of  axial  c.g.  dis¬ 
placement,  velocity,  and  acceler¬ 
ation  in  one  graph;  if  =  i 
curves  are  plotted;  if  =  0  or 

blank,  no  plot. 

Plot  curves  of  normal  c.g.  dis¬ 

placement,  velocity,  and  acceler¬ 
ation  in  one  graph;  if  =  1  curves 
are  plotted;  if  -  0  or  blank,  no 
plot . 

Plot  curves  of  normal  forces  at 

bourrelet  and  rotating  band;  if  = 
1,  plot  curves;  if  =  0  or  blank, 
no  plot. 

Plot  curve  of  normal  acceleration 
at  axial  point  F  of  projectile; 
if  =  1,  plot  curve;  if  =  0  or 

blank,  no  plot. 

Plot  curve  of  aerodynamic  drag 
force;  if  =  1,  plot  curve;  if  ==  0 
or  blank,  no  plot. 

Plot  curves  of  yaw  angle  and 

velocity;  if  =  1 ,  plot  curves;  if 
=  0  or  blank,  no  plot. 
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61-70 


NYAWA 


110 


Plot  curve  of  yaw  acceleration; 
if  =  1,  plot  curve;  if  -  0  or 

blank,  no  plot, 

71-80  NFGVSX  110  Plot  curves  versus  projectile 

travel;  if  =  0  or  blank,  curves 
are  plotted  versus  time  only;  if 
=  1,  curves  are  plotted  versus 

projectile  travel  in  addition  to 
those  plotted  versus  time. 

The  following  cards,  B  through  N,  are  omitted  if  no  figures  are  required - 
Card  B  Figure  scale  card 


List  MXMN ,  FGXMAXM,  FGTMAX,  FGXMAX 


Format  (110,  3F10.5) 


Columns  Symbol  Format 


Description 


1-10 


11-20 


21-30 


31-40 


MXMN 

110 

Determine  maximum  and  minimum 
values  of  graph  axis  scales;  if  = 
0  or  blank,  scale  values  will  be 
computed  in  the  program;  if  =  1  , 
values  are  entered  in  following 
cards  C-G. 

FGXMAXM 

F10.5 

Maximum  value  of  figure  displace¬ 
ment  axis  scale  in  metric  unit; 
if  left  blank,  it  will  be  com¬ 
puted  from  FGXMAX  value  in  the 
program. 

FGTMAX 

F10.5 

Maximum  value  of  figure  time  axis 
scale  in  ms;  if  =  0  or  blank  it 
will  be  set  to  equal  to  TMAX 
input  value  in  the  program, 

FGXMAX 

F10.5 

Maximum  value  of  figure  projec¬ 
tile  travel  axis  scale  in  English 
unit,  in;  leave  blank  if  FGXMAXM 
value  is  entered  or  if  FGXMAXM  is 
to  be  determined  in  the  program. 

If  MXMN  =  0  in  card  B,  the  following  cards  C  through  G  are  omitted. 


Card  C  Horizontal  axis  scale  and  unit  control  card 


List  FTMAX,  FTMIN ,  FXMAX,  FXMIN ,  UNITMS 

Format  (5F10.4) 
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Columns 


1-10 

11-20 

21-30 

31-40 

41-50 


Card  D 


List 


Format 


Columns 


1-10 


11-20 


21-30 


31-40 


Symbol  Format  Description 


FTMAX  F10.4 

FTMIN  F10.4 

FXMAX  F10.4 

FXMIN  F10.4 

UNITMS  FI 0.4 


Maximum  value  of  time  scale  of 
figure,  ms.  It  overrides  similar 
value  entered  in  card  B. 

Minimum  value  of  figure  time 
scale,  ms. 

Maximum  value  of  figure  travel 
axis  scale;  in  inches  if  UNITMS  = 
0  or  blank;  m  if  UNITMS  =  1. 

Minimum  value  of  figure  travel 
axis  scale;  in  inches  if  UNITMS  = 
0  or  blank;  m  if  UNITMS  =>  1. 

Determine  units  of  maximum  and 
minimum  values  of  axis  scale;  if 
=  0  or  blank,  English  units  are 
used;  if  =  1,  metric  units  are 

used. 


Axial  displacement,  velocity,  and  acceleration  scale  card 


FCGADMX,  FCGADMN,  FCGAVMX ,  FCGAVMN,  FCGAAMX ,  FCGAAMN 
(6F10.4) 

Symbol  Format  Description 

FCGADMX  F10.4  Maximum  value  of  projectile  axial 

displacement  axis  scale;  in 
inches  if  UNITMS  =  0  or  blank;  m 
if  UNITMS  =  1 . 


FCGADMN  F10.4 


FCGAVMX  FI 0.4 


FCGAVMN  FI 0.4 


Minimum  value  of  projectile  axial 
displacement  axis  scale;  in 
inches  if  UNITMS  *  0  or  blank;  m 
if  UNITMS  =1. 

Maximum  value  of  projectile  axial 
velocity  axis  scale;  in  ft/s  if 
UNITMS  =  0  or  blank;  m/s  if 
UNITMS  =  1. 

Minimum  value  of  projectile  axial 
velocity  axis  scale;  in  ft/s  if 
UNITMS  =  0  or  blank;  m/s  if 
UNITMS  =1. 
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41-50 

FCGAAMX 

F10.4 

Maximum  value  of  projectile  axial 
acceleration  axis  scale  in  g*s. 

51-60 

FCGAAMN 

F10.4 

Minimum  value  of  projectile  axial 
acceleration  axis  scale  In  g's. 

Card  E 

Normal  displacement, 

velocity,  and  acceleration  scale  card 

List 

FCGNDMX, 

FCGNDMN, 

FCGNVMX,  FCGNVMN,  FCGNAMX,  FCGNAMN 

Format 

(6F10.4) 

Columns 

Symbol 

Format 

Description 

1-10 

FCGNDMX 

F10.4 

Maximum  value  of  projectile  nor¬ 
mal  displacement  axis  scale;  in 
inches  if  UNITMS  =  0  or  blank;  mm 
if  UNITMS  =  1 . 

11-20 

FCGNDMN 

F10.4 

Minimum  value  of  projectile  nor¬ 
mal  displacement  axis  scale;  in 
inches  if  UNITMS  =  0  or  blank;  mm 
if  UNITMS  =  1. 

21-30 

FCGNVMX 

F10.4 

Maximum  value  of  projectile  nor¬ 
mal  velocity  axis  scale;  in  ft/s 
if  UNITMS  =  0  or  blank;  m/s  if 
UNITMS  =  1. 

31-40 

FCGNVMN 

F10.4 

Minimum  value  of  projectile  nor¬ 
mal  velocity  axis  scale;  in  ft/s 
if  UNITMS  =  0  or  blank;  m/s  if 
UNITMS  =  1. 

41-50 

FCGNAMX 

F10.4 

Maximum  value  of  projectile  nor¬ 
mal  acceleration  in  gfs. 

51-60 

FCGNAMN 

F10.4 

Minimum  value  of  projectile  nor¬ 
mal  acceleration  in  gfs. 

Card  F  Force  and  point  F  acceleration 


List  FLFMX,  FLFMN,  FAPFMX,  FAPFMN,  FAFMX,  FAFMN 

Format  (6F10.4) 

Columns  Symbol  Format  Description 

1-10  FLFMX  F10.4  Maximum  value  of  lateral  force 

axis  scale;  in  klb  if  UNITMS  =  0 

or  blank;  kN  if  UNITMS  «  1. 
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11-20 

FLFMN 

F10.4 

Minimum  value  of  lateral  force 
axis  scale;  in  klb  if  UNITMS  =  0 
or  blank;  kN  if  UNITMS  =  1. 

21-30 

FAPFMX 

F10.4 

Maximum  value  of  axial  point  F 
normal  acceleration  axis  scale  in 

g's. 

31-40 

FAPFMN 

F10.4 

Minimum  value  of  axial  point  F 
normal  acceleration  axis  scale  in 

g's. 

41-50 

FAFMX 

F10.4 

Maximum  value  of  air  resistance 
(aerodynamic  drag  force)  axis 
scale;  in  klb  if  UNITMS  =  0  or 
blank;  kN  if  UNITMS  =  1 . 

51-60 

FAFMN 

F10.4 

Minimum  value  of  air  resistance 
axis  scale;  in  klb  if  UNITMS  =  0 
or  blank;  kN  if  UNITMS  =  l . 

Card  G 

Yaw  angle,  velocity, 

and  acceleration 

List 

FYAWDMX, 

FYAWDMN, 

FYAWVMX,  FYAWVMN,  FYAWAMX,  FYAWAMN 

Format 

(6F10.4) 

Columns 

Symbol 

Format 

Description 

1-10 

FYAWDMX 

F10.4 

Maximum  value  of  yaw  angle  axis 
scale  in  milliradians. 

11-20 

FYAWDMN 

F10.4 

Minimum  value  of  yaw  angle  axis 
scale  in  milliradians. 

21-30 

FYAWVMX 

F10.4 

Maximum  value  of  yaw  velocity 
axis  scale  in  rad/s. 

31-40 

FYAWVMN 

F10.4 

Minimum  value  of  yaw  velocity 
axis  scale  in  rad/s 

41-50 

FYAWAMX 

F10.4 

Maximum  value  of  yaw  acceleration 
axis  scale  in  rad/s^. 

51-60 

FYAWAMN 

F10.4 

Minimum  value  of  yaw  acceleration 
axis  scale  in  rad/s^. 

Card  H 

Figure  subtitle  control  card 

List 

NT(I)  I  =  1,6 

25 


Format 


(6110) 


Columns 

1-10 


11-20 


21-30 


31-40 


41-50 


51-60 


Card  I 
If  NT(  1 ) 


Symbol  Format  Description 


NT(  1 )  110 


NT( 2)  110 


NT( 3)  110 


NT( 4)  110 


NT( 5)  110 


NT(  6)  no 


Control  of  subtitle  1;  if  =  1, 

following  subtitle  1  card  (card 

I)  will  be  read;  0  or  blank, 
subtitle  1  card  is  omitted  and 
title  of  analysis  (TITLEM)  will 
be  used  as  subtitle  1  in  output 
figures. 

Control  of  subtitle  2;  if  *  1  , 

following  subtitle  2  card  (card 

J)  will  be  read;  0  or  blank, 
subtitle  2  card  is  omitted  and 
subtitle  of  analysis  (TITLES) 
will  be  used  as  subtitle  2  in 
output  figures. 

Control  of  subtitle  3;  if  =  1  , 

following  subtitle  3  card  (card 

K)  will  be  read;  0  or  blank, 
subtitle  3  card  is  omitted  and 
program  will  print  c.g.  eccentri¬ 
city  data  as  subtitle  3. 

Control  of  subtitle  4;  if  =  1 , 

following  subtitle  4  card  (card 

L)  will  be  read;  0  or  blank, 
subtitle  4  card  is  omitted  and 
program  will  print  bore  and  bour- 
relet  diameter  data  as  subtitle 
4. 

Control  of  subtitle  5;  if  =  1 , 

following  subtitle  5  card  (card 

M)  will  be  read;  0  or  blank, 
subtitle  5  card  is  omitted  and 
program  will  print  friction  coef¬ 
ficient  data  as  subtitle  5. 

Control  of  subtitle  6;  if  =  1, 

following  subtitle  6  card  (card 

N)  will  be  read;  0  or  blank, 
subtitle  6  card  is  omitted  and 
program  will  print  spring  con¬ 
stant  data  as  subtitle  6. 


Subtitle  1  card 


=  0  or  blank,  this  card  is  omitted. 
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List 


Format 


Columns 


TITLR(J,1)  J=1 ,5 
( 5A10) 

Symbol  Format 


Description 


1-50 

TITLE( J , 1 )  5A10 

Subtitle  1  to  be 
put  figures. 

printed 

in 

out- 

Card  J 

Subtitle  2  card 

If  NT( 2) 

=  0  or  blank,  this  card 

is  omitted. 

List 

TITLE( J,2)  J-1,5 

Format 

( 5A10) 

Columns 

Symbol  Format 

Description 

1-50 

TITLE( J ,2)  5A10 

Subtitle  2  to  be 
put  figures. 

printed 

in 

out- 

Card  K 

Subtitle  3  card 

If  NT( 3) 

=  0  or  blank,  this  card 

is  omitted. 

List 

TITLE( J,3)  J-1,5 

Format 

( 5A10) 

Columns 

Symbol  Format 

Description 

1-50 

TITLE( J ,3)  5A10 

Subtitle  3  to  be 
put  figures. 

printed 

in 

out- 

Card  L  Subtitle  4  card 

If  NT(4) 

=  0  or  blank,  omit  this 

card . 

List 

TITLE( J ,4)  J=1 ,5 

Format 

(5A10) 

Columns 

Symbol  Format 

Description 

1-50 

TITLE( J,4)  5A10 

Subtitle  4  to  be 
put  figures. 

printed 

in 

out- 

Card  M  Subtitle  5  card 

If  NT(5)  -  0  or  blank,  omit  this  card. 


27 


List 

TITLE( J,5) 

J=1 ,5 

Format 

(5A10) 

Columns 

Symbol 

Format 

Description 

1-50 

TITLE( J ,5) 

5A10 

Subtitle  5  to  be  printed  in  out¬ 
put  figures. 

Card  N 

Subtitle  6  card 

If  NT( 6) 

=  0  or  blank, 

omit  this 

card. 

List 

TITLE( J ,6) 

J-1,5 

Format 

( 5A10) 

Columns 

Symbol 

Format 

Description 

1-50  TITLE(J,6)  5A10  Subtitle  6  to  be  printed  in  out¬ 

put  figures. 

9.  Firing  pressure  card  set  (2-26  cards) 

Card  A  General  information  card 


List  NOFP,  NCHAMP,  CWT,  PCGX,  UNITFP 


Format  (2110,  4F10.5) 


Columns  Symbol  Format 

1-10  NOFP  110 

11-20  NCHAMP  110 

21-30  CWT  F10.5 

31-40  PCGX  F10.5 

41-50  UNITFP  F10.5 


Description 

Number  of  pressure-time  data 
pairs . 

Type  of  firing  pressure;  enter  1 
for  chamber  pressure;  0  or  blank 
for  base  pressure. 

Weight  of  propellant  charge;  in 
lb  if  UNITFP  =  0  or  blank;  kg  if 
UNITFP  -  1. 

Base  pressure  center  to  c.g. 
distance;  in  inches  if  UNITFP  =  0 
or  blank;  cm  if  UNITFP  =  1.0; 
default  value  =  CGX. 

Unit  of  pressure  input;  if  =  0  or 
blank,  pressure  in  psi;  if  =  1.0, 
pressure  in  kPa. 
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Card  B  Time-pressure  cards 
List  FPT(I),  FP(I),  1=1,  NOFP 

(4(F10.5,  F10 . 2) ) 

Symbol  Format 

FPT(I) 

FP(I) 


Format 


Columns 


1-10 

11-20 


F10.5 

F10.2 


Description 
Time  in  ms 


Pressure  at  time  FPT(l);  in  psi 
if  UNITFP  =  0  or  blank;  kPa  if 
UNITFP  =1.0. 


Enter  successive  time-pressure  data  pairs  in  successive  columns:  21-30,  31-40, 
41-50,  51-60,  61-70,  and  71-80.  Repeat  card  B  procedure  until  all  data  are  en¬ 
tered. 


10.  Worn  tube  card  set  (1-26  cards) 


Card  A  General  information  card 


List 

NOTW,  UNITTW 

Format 

(110,  F10.5) 

Columns 

Symbol  Format 

Description 

1-10 

NOTW  110 

Number  of  tube  wear  data 

pairs . 

11-20 

UNITTW  F10.5 

Unit  of  tube  wear  data; 

if  =  0  or 

blank,  data  in  English  units;  if 
=  1.0,  data  in  metric  units. 


If  NOTW  =  0,  following  card  B  is  omitted. 


Card  B  Distance-wear  cards 


List  TWMX(I),  TWDC(I),  1=1,  NOTW 

Format  (8F10.5) 


Columns 

Symbol 

Format 

1-10 

TWMX(I) 

F10.5 

11-20 

TWDC(I) 

F10.5 

Description 

Distance  from  muzzle  to  tube  wear 
location,  in  inches  if  UNITTW  =  0 
or  blank;  cm  if  UNITTW  =  1.0. 

Change  of  bore  diameter;  in  mils 
of  inch  if  UNITTW  =  0  or  blank; 
mm  if  UNITTW  =  1.0. 
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Enter  successive  data  pairs  in  successive  columns:  21-30,  31-40,  41-50,  51-60 
61-70,  and  71-80.  Repeat  card  B  procedure  until  all  data  are  entered. 


11.  Bourrelet  elasticity  card  set  (1-26  cards) 

Card  A  General  information  card 


List  NOEBK,  BRTK,  BRTE ,  UNITBK 


Format  (HO,  3F10.4) 


Columns 

1-10 

11-20 


Symbol  Format 

NOEBK  HO 

BRTK  F10.4 


Description 


Number  of  experimental  bourrelet 
deflection-load  data  pairs. 

Average  or  linear  bourrelet 
spring  constant,  used  when  exper¬ 
imental  deflection-load  data  are 
not  available;  in  106  lb/in  if 
UNITBK  -  0  or  blank;  106  N/cm  if 
UNITBK  =  1.0. 


21-30 


31-40 


BRTE  F10.4  Young’s  modulus  of  bourrelet 

material ,  used  to  compute  approx¬ 
imate  spring  constant  when  test 
data  and  average  spring  constant 
are  not  available;  in  106  psi  if 
UNITBK  =  0  or  blank;  10^  kPa,  if 
UNITBK  =1.0. 

UNITBK  F10.4  Unit  of  bourrelet  elasticity 

data;  if  =  0  or  blank,  data  in 
English  units;  if  =  1.0,  data  in 
metric  units. 


If  NOEBK  -  0,  omit  following  card  B. 


Card  B  Deflect  ion— load  cards 


List 


BDEF(I),  BLOAD(I),  1=1,  NOEBK 


Format  (4(F10.5,  F10.2)) 


Columns 


Symbol  Format  Description 


BDEF(I)  F10.5  Diametrical  deflection  of  bourre— 

let;  in  inches  if  UNITBK  =  0  or 
blank;  cm  if  UNITBK  =  1.0. 
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11-20 


BLOAD(I)  F10.2 


Load  applied  at  bourrelet  deflec¬ 
tion  test  corresponding  to  de¬ 
flection  BDEF(I);  in  lb  if  UNITBK 
=  0  or  blank;  N  if  UNITBK  =1.0. 


Enter  successive  deflection-load  data  pairs  in  successive  columns:  21-30,  31-40, 
41-50,  51-60,  61-70,  and  71-80.  Repeat  card  B  procedure  until  all  data  are  en¬ 
tered  • 


12.  Rotating  band  elasticity  card  set  (1-51  cards) 


Card  A  General  information  card 


List  NOERBK,  NOERBMK,  RBK,  RBMK,  RTBAND ,  UNITRBK 

Format  (2110,  4F10.4) 

Columns  Symbol  Format  Description 


1-10  NOERBK  110 


Number  of  experimental  rotating 
band  deflection  test  data  pairs. 


11-20 


21-30 


NOERBMK  110 


RBK  F10.4 


Number  of  experimental  transverse 
rotation  at  rotating  band  test 
data  pairs. 


Average  or  linear  spring 
of  rotating  band;  in  10° 
UNITRBK  =  0  or  blank;  106 


constant 
lb/in  if 
N/cm  if 


UNITRBK  =  1.0. 


31-40 


41-50 


51-60 


RBMK  F10.4 


RTBAND  F10.4 


UNITRBK  FI 0.4 


Average  or  linear  angular  spring 
constant  of  rotating  band;  in  10° 
in-lb/deg  if  UNITRBK  =  0  or 

blank;  106  cm-N/rad  if  UNITRBK  = 

1.0, 


Identification  of  rotating  band; 
=1.0  when  rotating  band  is  used; 
=  0  or  blank,  when  obturator  is 
used . 

Unit  of  spring  constant  data;  if 
=  0  or  blank,  English  units  are 
used;  if  =  1.0,  metric  units  are 
used . 


If  NOERBK  =  0,  omit  following  card  B. 

Card  B  Deflection-load  cards 


List  RBDEF(I),  RBLOAD(I),  1=1,  NOERBK 
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Format 

(4(F10.5, 

F10 . 2) ) 

Columns 

Symbol 

Format 

Description 

1-10 

RBDEF(I) 

F10.5 

Diametrical  deflection  of  rotat¬ 
ing  band  in  inches  if  UNITRBK  =  0 
or  blank;  cm  if  UNITRBK  =  1.0. 

11-20 

RBLOAD(I) 

F10.2 

Load  applied  at  rotating  band 
deflection  test  corresponding  to 
deflection  RBDEF(I) ;  in  lb  if 

UNITRBK  =  0  or  blank;  N  If 

UNITRBK  =  1.0. 

Enter  successive  deflection-load  data  pairs  in  successive  columns:  21-30,  31-40, 
41-50,  51-60,  61-70,  and  71-80.  Repeat  card  B  procedure  until  all  data  are  en¬ 
tered.  If  NOERBMK  =  0  or  blank,  omit  following  card  C. 

Card  C  Angular  deflection— moment  cards 


list  RBANG(I),  RBM(I),  1=1,  NOERBMK 


Format  (4(F10.5,  F10.2)) 

Columns  Symbol  Format 

1-10  RBANG(I)  F10.5 

11-20  RBM(I)  F10.2 


Description 


Angular  deflection  or  transverse 
rotation  of  projectile  at  rotat¬ 
ing  band;  in  deg  if  UNITRBK  =  0 
or  blank;  rad  if  UNITRBK  =  1.0. 

Moment  applied  at  rotating  band 
angular  deflection  test  corre¬ 
sponding  to  angular  deflection 
RBANG(I);  in  in-lb  if  UNITRBK  =  0 
or  blank;  cm-N,  if  UNITRBK  =  1.0. 


Enter  successive  angular  deflection-moment  data  pairs  in  successive  columns:  21- 
30,  31-40,  41-50,  51-60,  61-70,  and  71-80.  Repeat  card  C  procedure  until  all 
data  are  entered. 


13.  Friction  card  set  (1-27  cards) 

Card  A  Friction  coefficient  and  general  information  card 

List  NORBPR,  RBPRAV ,  RBMU ,  BRTMU ,  RBPW,  RBPRMR,  RBNU 

UNITRBP 

Format  (110,  F10.2,  6F10.5) 
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Columns 

Symbol 

Format 

Description 

1-10 

NORBPR 

no 

Number  of  gun  tube-rotating 
contact  pressure  data  pairs. 

band 

11-20 

RBPRAV 

F10.2 

Average  gun  tube-rotating 
contact  pressure;  in  psi 

band 

if 

UNITRBP  =  0  or  blank;  kPa  if 
UNITRBP  =  1.0;  not  entered  when 

NORBPR  is  not  zero  or  when  pres¬ 
sure  will  be  computed  in  the 
program. 


21-30 

RBMU 

F10.5 

Friction  coefficient  at  rotating 
band . 

31-40 

BRTMU 

F10.5 

Friction  coefficient  at  bourre- 
let . 

41-50 

RBPW 

F10.5 

Effective  width  of  rotating  band 
for  gun  tube-rotating  band  con¬ 
tact  pressure  computation  in  the 
program;  in  inches  if  UNITRBP  =  0 
or  blank;  cm  if  UNITRBP  *  1.0; 

default  value  =  RBW. 

51-60 

RBPRMR 

F10.5 

Known  or  estimated  ratio  of  gun 
tube-rotating  band  contact  pres¬ 
sure  at  muzzle  to  that  at  breech; 
default  value  =  0.5. 

61-70 

RBNU 

F10.5 

Poisson’s  ratio  of  rotating  band 
material • 

71-80 

UNITRBP 

F10.5 

Unit  of  gun  tube-rotating  band 
contact  pressure  input  data;  =  0 
or  blank  if  English  units  are 
used;  =  1.0  if  metric  units  are 

used. 

If  RBPRAV  is  not  zero,  omit  following  cards  B  and  C.  If  RBNU  is  zero,  omit  fol¬ 
lowing  card  B. 


Card  B  Rotating  band  card 


List  RBOD ,  RBID,  RBE ,  PROJRBD 

Format  (4F10.5) 
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Columns 

Symbol 

Format 

1-10 

RBOD 

F10.5 

11-20 

RBID 

F10.5 

21-30  RBE  F10.5 


31-40  PROJRBD  F10.5 


Description 

Initial  outer  diameter  of  rotat¬ 
ing  band;  in  inches  if  UNITRBP  = 
0  or  blank;  cm  if  UNITRBP  =  1.0 

Initial  inner  diameter  of  rotat¬ 
ing  band  if  it  is  different  from 
projectile  diameter  at  band  seat, 
or  initial  inner  diameter  of 
obturator  before  it  is  fitted  to 
projectile;  in  inches  if  UNITRBP 
-  0  or  blank;  cm  if  UNITRBP  = 
1.0;  default  value  —  band  seat 
diameter  of  projectile  (PROJRBD) 

Young’s  modulus  of  rotating  band 
material;  in  10°  psi  if  UNITRBP  = 
0  or  blank;  106  kPa  if  UNITRBP  = 

1.0. 

Diameter  of  band  seat  of  projec¬ 
tile;  in  inches  if  UNITRBP  =  0  or 
blank;  cm  if  UNITRBP  =  1.0. 


If  NORBPR  is  zero,  omit  folliwng  card  C. 


Card  C  Distance-pressure  card 


List 


RBPMX(I),  RBPR(I),  1=1,  NORBPR 


Format  (4(F10.5,  F10.2)) 


Columns 

Symbol 

Format 

1-10 

RBPMX(I) 

F10.5 

11-20 

RBPR(I) 

F10.2 

Description 

Distance  of  gun  tube-rotating 
band  pressure  location  from  muz¬ 
zle;  in  inches  if  UNITRBP  =  0  or 
blank;  cm  if  UNITRBP  =  1.0. 

Gun  tube-rotating  band  pressure; 
in  psi  if  UNITRBP  =  0  or  blank; 
kPa  if  UNITRBP  =1.0. 


Enter  successive 
40,  41-50,  51-60 
entered . 


distance-pressure  data  pairs  in  successive  columns:  21-30, 
61-70,  and  71-80.  Repeat  card  C  procedure  until  all  data 
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are 


14.  Air  resistance  card  set  (1-77  cards) 

Card  A  General  information  and  drag  card 


34 


List 


NOCD,  CDAV,  RHO,  AC,  DTALAM,  UNTTAF 


Format 

(110,  5F10 

.5) 

Columns 

Symbol 

Format 

1-10 

NOCD 

110 

11-20 

CDAV 

F10.5 

21-30 

RHO 

F10.5 

31-40  AC  F10.5 


41-50  DTALAM  F10.5 


Description 


Number  of  yaw  angle-aerodynamic 
drag  coefficient  data  pairs. 


Average  or  constant  drag  coeffi¬ 
cient  . 


Air  mass  density;  in  lbm/in'*  if 
UNITAF  =  0  or  blank;  kg/cm^  if 

UNITAF  =  1.0;  default  value  = 

standard  density  at  sea  level 
(0.002378  slug/ft^)  is  used  in 
the  program. 


Location  of  aerodynamic  center  or 
axial  distance  from  projectile 
c.g.  section  to  axial  point  where 
aerodynamic  pressure  center  in¬ 
tersects;  in  inches  if  UNITAF  =  0 
or  blank;  cm  if  UNITAF  *1.0. 


Control  number  of  aerodynamic 
lift  and  moment;  if  *  0  or  blank, 
no  data  of  lift  and  moment  are 
used  in  the  computation  or  lift 
and  moment  are  ignored;  if  =  1.0, 
data  are  entered  in  following 
cards  C  and  D. 


51-60  UNITAF  FI 0.5 


Unit  of  aerodynamic  yaw- 
coefficient  data;  data  in  English 
units  if  UNITAF  =  0  or  blank;  in 
metric  units  if  UNITAF  -  1.0. 


If  NOCD  =  0,  omit  following  card  B. 


Card  B  Drag  coefficient  card 

List  CDANG(I),  CD(I),  1*1,  NOCD 

Format  (8F10.5) 

Columns  Symbol  Format  Description 

1-10  CDANG(I)  F10.5  Yaw  angle  or  angle  of  attack;  in 

deg  if  UNITAF  =  0  or  blank;  rad 

if  UNITAF  =1.0. 
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11-20 


CD  ( I ) 


F10.5 


Aerodynamic  drag 
corresponding  to 
CDANG(I). 


coefficient 
yaw  angle 


Enter  successive  angle-coefficient  data  pairs  in  successive  columns:  21-30,  31- 
40,  41-50,  51-60,  61-70,  and  71-80.  Repeat  card  B  procedure  until  all  data  are 
entered.  If  DTALAM  =  0  or  blank,  omit  following  cards  C-E. 

Card  C  Lift  and  moment  information  card 


List 

NOCL, 

NOCM,  CLAV, 

CMAV,  CHORD 

Format 

(2110, 

3F10.5) 

Columns 

Symbol 

Format 

Description 

1-10 

NOCL 

110 

Number  of  aerodynamic  lift  coef¬ 
ficient  data  pairs. 

11-20 

NOCM 

110 

Number  of  aerodynamic  moment 

coefficient  data  pairs. 

21-30 

CLAV 

F10.5 

Average  or  constant  aerodynamic 
lift  coefficient;  entered  when 
NOCL  =  0. 

31-40 

CMAV 

F10.5 

Average  or  constant  aerodynamic 
moment  coefficient;  entered  when 
NOCM  «  0. 

41-50 

CHORD 

F10.5 

Length  of  equivalent  airfoil 

chord  for  aerodynamic  moment 
computation;  in  inches  if  UNITAF 
=  0  or  blank;  cm  if  UNITAF  =  1.0; 
default  value  =  PROJL  (length  of 
projectile) . 

If  N0CL  =  0,  omit  following  card  D. 

Card  D  Lift  coefficient  card 

List  CLANG(I),  CL( I) ,  1=1,  NOCL 

Format  (8F10.5) 

Columns  Symbol  Format  Description 

l~lt>  CLANG(I)  F10.5  Yaw  angle  or  angle  of  attack;  in 

deg  if  UNITAF  =  0  or  blank;  rad 

if  UNITAF  =1.0. 
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11-20  CL(I)  F10.5  Aerodynamic  lift  coefficient 

corresponding  to  CLANG(I). 

Enter  successive  angle-coefficient  data  pairs  in  successive  columns:  21-30,  31- 
40,  41-50,  51-60,  61-70,  and  71-80.  Repeat  card  D  procedure  until  all  data  are 
entered. 


Card  E  Moment  coefficient  card 


If  NOCM  =  0,  omit  this  card. 


List  CMANG(I),  CM(I),  1=1,  NOCM 


Format  (8F10.5) 


Columns 

Symbol 

Format 

1-10 

CMANG(I) 

F10.5 

11-20 

CM(I) 

F10.5 

Description 

Yaw  angle  or  angle  of  attack;  in 
deg  if  UNITAF  =  0  or  blank;  rad 
if  UNITAF  =  1.0. 

Aerodynamic  moment  coefficient 
corresponding  to  CMANG(I). 


Enter  successive  angle-coefficient  data  pairs  in  successive  columns:  21-30,  31- 
40,  41-50,  51-60,  61-70,  and  71-80.  Repeat  card  E  procedure  until  all  data  are 
entered. 


15.  End  card  (1  card) 

An  end-of-record  (EOR)  job  control  (7-8-9)  card  which  consists  of  a 
7-8-9  multipunch  in  column  1  is  used  as  the  end  card  of  the  input  data  deck. 

To  facilitate  the  preparation  of  the  input  data  cards  and  for  easy 
reference,  the  above-mentioned  cards  are  summarized  in  appendix  A. 


OPERATION 


The  program  is  operated  with  job  control  and  data  cards.  It  has  been  oper¬ 
ated  on  a  CDC  6000  series  computer  with  an  in-line  Nicolet  Zeta  digital  drum 
plotter  or  similar  one  for  curve  plotting.  It  may  be  adapted  to  be  operated  at  a 
Tektronix  4014  or  similar  graphics  terminal  to  get  an  ’’instant"  copy  of  a  fig¬ 
ure.  The  program  requires  approximately  120000g  words  of  core.  Detailed  infor¬ 
mation  of  job  deck  setup  for  the  computer  run  may  be  obtained  from  the  operator 
of  the  user’s  computer  terminal  and  is  not  discussed  here. 

It  should  be  noted  that  the  accuracy  of  computation  will  depend  on  the  input 
values  of  different  data,  especially  the  spring  constants  and  other  force  and 
moment  coefficients.  Experimental  data  are  preferred  whenever  they  are  availa- 
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ble.  Furthermore,  the  time  step  size  affects  the  processing  and  the  accuracy  of 
numerical  solution.  Therefore,  it  is  important  to  use  very  small  time  step  sizes 
(microseconds  or  less)  to  obtain  good  results.  The  F  conversion  formats  used  in 
the  program  are  suitable  for  general  projectile  computations.  To  accommodate  for 
larger  or  smaller  values  they  should  be  changed  to  E  conversion  formats  if  neces- 

S  3.  Tl'Y  m 


EXAMPLE  OF  COMPUTATION 


An  example  of  using  the  program  is  presented  here.  The  problem  is  to  anal- 
yz  a  7  kg  projectile  launched  in  a  120-mm  smooth  gun  tube.  Most  of  the  input 
data,  such  as  the  dimensions  and  physical  properties  of  the  projectile  and  gun 
tube,  are  obtained  from  design  drawings,  design  error  budget,  or  data  of  similar 
projectiles.  Other  data  are  determined  from  experimental  results,  approximate 
computations,  or  engineering  handbooks. 

tile  kfilSllle  C*-f;  eccentricity  has  a  magnitude  of  0.254  mm  and  the  projec- 
tile  !s  initially  positioned  so  that  the  c.g.  is  at  the  six  o'clock  position  and 
the  bourrelet  contacts  the  gun  tube.  Thus,  the  c.g.  is  under  the  projectile 
geometrical  axis.  It  is,  therefore,  taken  as  -0.254  mm  for  the  input  value.  The 
six  o  clock  position  indicates  that  the  X-Y  plane  is  vertical  and,  therefore  8 

s  zero.  The  gun  tube  is  held  in  a  horizontal  position,  therefore,  a  is  also 
zero  • 


of  tl’e  bou"elet  ls  obtained  from  the  compression  tests 

,  ,th?  Projectile,  and  Che  average  value  of  the  results  Is  used.  The  bourrelet 
oa  deflection  curve  data  may  be  used  instead  of  the  average  value  if  more  accu- 
rate  results  are  desired.  The  angular  spring  constant  of  the  obturator  is  taken 
zero  since  the  yaw  angle  is  small  and  no  test  data  are  available.  The  fric- 
onal  coefficients  and  the  in-bore  aerodynamic  coefficients  are  also  not  avail¬ 
able  at  present.  The  lift  and  moment  coefficients  of  air  resistance  are  taken  as 
zero  since  they  are  small.  The  drag  coefficient  is  assumed  to  be  unity  in  order 
o  ge  some  idea  of  the  magnitude  of  the  in-bore  air  resistance.  The  frictional 
coefficients  at  the  bourrelet  and  the  obturator  are  assumed  to  be  0.15  and  0.04 
based  on  the  similar  data  of  steel  to  steel  and  teflon  to  steel  (ref  8). 

The  base  firing  pressure-time  data  are  determined  from  an  interior  ballis¬ 
tics  computation  as  shown  in  figure  23  and  in  table  3  of  appendix  B. 

c,  „The(-dKtailedJ  iJn,pUt  dimensions  and  physical  properties  of  the  projectile  and 
the  gun  tube  are  indicated  in  table  1  of  appendix  B.  J 

BAT  TOT?6"  tThle  computation  is  "Example  of  Computation  of  Program 

uLO  2.  he  subtitle  is  blank.  English  units  are  used  in  the  input  data  pre¬ 
paration.  Therefore,  NDAMS  and  NDARAD  are  both  zero.  Other  unit  control  numbers 
except  NPOPDTA  (=20)  are  also  zero.  Consequently,  the  default  values  In 
program  are  used.  Then  the  input  data  tables  are  printed  out;  the  intermediate 
computation  check  values  are  printed  out  every  15th  value  recorded  on  the  magne- 
c  tape;  the  output  table  of  maximum,  minimum,  and  muzzle  values  is  printed;  and 
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every  20th  value  of  the  computed  values  recorded  on  the  magnetic  tape  is  printed 
out. 

The  maximum  time  of  computation  is  4.8  ms  and  the  step  size  is  0.001  ms. 
The  constant  step  size  computation  is  required  so  TSLAST  is  zero. 

All  output  figures  plotted  against  both  time  and  travel  are  desired,  conse¬ 
quently,  all  figure  plot  controls  are  1.  The  maximum  and  minimum  scales  of  out¬ 
put  figures  are  to  be  determined  by  the  program  and,  therefore,  are  all  zero. 

The  title  of  this  computation  will  be  used  as  one  of  the  subtitles  of  the 
output  figures  and  no  special  subtitles  are  desired.  Therefore,  all  NT(I)  are 
zero  and  no  subtitle  cards  are  needed. 

Other  input  data  are  determined  according  to  instructions  stated  in  the 
input  data  coding.  The  input  data  card  arrangement  is  easily  understood  by  com¬ 
paring  it  (table  1  of  appendix  B)  with  the  input  data  tables  (tables  2  through  5 
of  appendix  B). 

With  the  input  data  determined,  the  input  data  cards  are  punched  according 
to  the  coding  formats  mentioned  before.  These  data  cards,  together  with  the 
program  and  job  control  cards,  form  the  job  card  deck  which  is  placed  in  the  card 
reader  of  a  CDC  6000  series  computer  terminal  to  be  read  and  executed  by  the 
computer.  The  printed  output  and  plotted  figures  are  shown  in  appendix  B. 
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Table  1*  List  of  labeled  common  blocks 


Common  block 

CMXMN 

CMXMNA 

CMXMNT 

CMXMNX 

DIM 

DTDI 

DTDIMA 

DTDIMB 

DTDIMC 

DTDIMD 

DTDIOM 

DTDIOMF 

DTM 

FGCVA 

FGCVB 

FGCVC 

FGCVE 

FGCVF 

FGDIA 

FGDIB 

FGDIVM 

OPTBMA 

OPTBMB 


Used  in 


BALLOT2  DIPTBL  PATAINP  FIGCTRL  FGOV3  FGPA3  FGMN3  FGPRTSB 

X  X 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


X 

X 

X 


X 

X 

X 

X 

X 

X 

X 

X 


X 

X 

X 

X 

X 

X 

X 

X 


X 

X 

X 


X 

X 


X 

X 

X 

X 


X 

X 

X 


X 

X 


X 

X 

X 


OPTBL 

X 

X 

X 

X 


X 

X 


X 

X 
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Table  2. 


Definitions  of  symbols  used  in  main 


program 


Symbol  Description 


A 

Area  of  gun  tube  bore  -  in^ 

AC 

Axial  distance  of  projectile  c.g.  to  aerodynamic  center, 
used  for  air  resistance  computation  -  in 

ADD;  FDD 

Normal  acceleration  at  projectile  axial  point  A  F  - 
in/s  * 

ADDM,  ADDN; 

FDDM,  FDDN 

Maximum,  minimum  value  of  ADD; FDD 

ADDMT ,  ADDNT ; 
FDDMT,  FDDNT 

Time  when  ADD; FDD  is  maximum,  minimum  -  ms 

ADDMX ,  ADDNX ; 
FDDMX,  FDDNX 

Travel  of  projectile  when  ADD; FDD  is  maximum,  minimum  - 
in 

ALPHA 

Gun  tube  elevation  angle  -  rad 

AM 

Total  aerodynamic  moment  about  projectile  c.g.  -  in-lb 

APFXM 

Distance  of  projectile  axial  point  F  from  projectile  nose 
in  metric  unit  -  cm 

AR 

Ratio  of  specific  heats  of  propellant  gases 

BDEF 

Deflection  of  bourrelet  compression  test  -  in 

BETA 

Angle  between  vertical  and  x-y  plane  -  rad 

BK 

Average  spring  constant  of  bourrelet  -  lb/in 

BLOAD 

Load  applied  at  bourrelet  compression  test  -  lb 

BM 

Total  moment  of  all  forces  at  bourrelet  -  in-lb 

BMU ,  BRTMU 

Friction  coefficient  at  bourrelet 

BN 

Normal  force  at  bourrelet  -  lb 

BNM,  BNN 

Maximum,  minimum  of  BN 

BNMT,  BNNT; 

BNMX,  BNNX 

Time;  travel  of  projectile  when  BN  is  maximum,  minimum  - 
ms 

BNRM 

BN  expressed  in  newtons 
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Table  2.  (cont) 


Symbol  Description 


BORED 

Gun  tube  bore  diameter  -  in 

BOREDM 

Bore  diameter  expressed  in  cm 

BP 

Base  pressure  used  in  computation  -  psi 

BPSLP 

Slope  of  base  pressure-time  curve 

BRTE 

Young* s  modulus  of  bourrelet  -  10^  psi 

BRTK 

Average  spring  constant  of  bourrelet  -  106  lb/in 

BRTKM 

BRTK  expressed  in  metric  unit  -  10^  N/cm 

BRTOD 

Bourrelet  outside  diameter  -  in 

BRTODM 

BRTOD  expressed  in  metric  unit  -  cm 

BRTTK 

Thickness  of  projectile  wall  at  bourrelet  -  in 

BRTW 

Length  of  bourrelet  portion  having  diameter  =  BRTOD  -  in 

BRTX 

Distance  of  bourrlet  from  projectile  base  plane  -  in 

CD,  CL,  CM 

Aerodynamic  coefficient  of  drag,  lift,  moment  for  projec¬ 
tile  in-bore  motion 

CDANG ,  CLANG ,  CMANG 

Yaw  angle  for  computing  in-bore  aerodynamic  coefficient 
CD,  CL,  CM  -  rad 

CDAV,  CLAV,  CMAV 

Average  CD,  CL,  CM 

CGEL 

c.g.  eccentricity  expressed  in  distance  -  in 

CGELM 

CGEL  expressed  in  metric  unit  -  cm 

CGEMT 

Projectile  c.g.  eccentricity  expressed  in  unbalanced 
moment  -  in-oz 

CGEMTM 

CGEMT  expressed  in  metric  unit  -  cm-kg 

CGX 

Distance  of  projectile  c.g.  to  base  plane  -  in 

CGXM 

CGX  expressed  in  metric  unit  -  cm 

CHORD 

Reference  length  for  in— bore  aerodynamic  moment  computa¬ 
tion  -  in 
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Table  2.  (cont) 


Symbol  Description 


CTHETA 

Cosine  of  angle  theta 

CW 

Ratio  fo  charge  weight  to  projectile  weight 

CWT 

Charge  weight  -  lb 

Table  2.  (cont) 

DABSGN 

Sign  of  bourrelet  deflection 

DAOSGN 

Sign  of  rotating  band  deflection 

DELTAB ,  DELTAO 

Maximum  normal  deflection  at  bourrelet,  rotating  band  - 
in 

DELTAT 

Time  step  size  -  sec 

FB 

Frictional  force  at  bourrelet  -  lb 

FD,  FL 

Aerodynamic  drag,  lift  force  -  lb 

FLM,  FLN 

Maximum,  minimum  FL  -  lb 

FLMT,  FLNT 

Time  when  FL  is  maximum,  minimum  -  ms 

FLMX,  FLNX 

Travel  of  projectile  when  FL  is  maximim,  minimum  -  in 

FOl 

Frictional  force  at  rotating  band  due  to  uniform  gun 
tube  rotating  band  contact  pressure  —  lb 

F01CF 

Product  of  tube— band  contact  area  and  friction  coeffi¬ 
cient 

F02 

Friction  force  due  to  normal  force  at  rotating  band  -  lb 

FP 

Base  pressure  -  psi 

FPEA 

Effective  base  pressure  acting  area  -  in^ 

FPT 

Time  value  of  firing  pressure-time  data  -  ms 

FX,  FY 

Total  force  acting  at  projectile  c.g.  in  X,  Y  direction  - 
lb 

GX,  GY 

X,  Y  component  of  gravitational  force  at  projectile  c.g. 

—  lb 

H 

Distance  of  bourrelet  to  c.g.  section  -  in 
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Table  2.  (cont) 


Symbol 

HH 

HXL 

HXLM 

NBK,  NBRTKC 
NC 

NCCK 

NCO,  NCL,  NCM 

NCHAMP 

NFP,  NOK,  NOMK, 

NIDTBL,  NOOPDTA, 
NOPTBL,  NPOPDTA 

NOCD  ,  NOCL,  NOCM 
NOEBK,  NOERBK, 
NOERBMK,  NOFP, 
NORBPR,  NOTW 

NOCP 

NR 

NRD,  NWD 
NT 

NTBL 

NWRB,  NWRO 


Description 


Projectile  c.g.  distance  to  muzzle  -  in 

Distance  between  bourrelet  and  rotating  band  mid-section 
-  in 

HXL  expressed  in  metric  unit  -  cm 

Control  number  in  bourrelet  spring  constant  computation 
Total  number  of  computations 

Control  number  in  printing  out  some  computed  results 
during  computation 

Control  number  in  aerodynamic  coefficient  computation  for 
drag,  lift,  moment 

Control  number  in  conversion  of  chamber  to  base  pressure 

NOP  Control  number  in  base  pressure,  gun  tube-rotating  band 
contact  pressure,  angular  spring  constant,  given  tube- 
band  contact  pressure  computation 

Control  number  in  printing  out  input  data  table,  record¬ 
ing  output  data,  printing  out  output  data  table,  detailed 
output  data 

Number  of  input  data  pairs  for  input  data  of  aerodynamic 
drag,  lift,  moment,  bourrelet  spring  constant,  rotating 
band  spring  constant,  angular  spring  constant,  firing 
pressure,  gun  tube-rotating  band  contact  pressure,  gun 
tube  wear 

Estimated  total  number  of  computations 
Number  of  computed  results  recorded  on  tape  1 
Control  number  in  recording  computed  data 
Control  number  in  reading  figure  subtitle 
Number  of  tables  printed 

Control  number  in  computing  tube  wear  at  bourrelet,  ro¬ 
tating  band 
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Table  2.  (cont) 


Symbol  Description 


OK,  OMK 

Spring  constant  lb/in,  angular  spring  constant  of  ro¬ 

tating  band  -  in-lb/rad 

OMU,  RBMU 

Friction  coefficient  of  rotating  band 

0M1,  OM2 ,  0M3 

Moment  about  projectile  c.g.  from  rotating  band  friction 
force  ( FO 1 ) ,  normal  force  (ON),  and  friction  force  (F02) 
transverse  rotation  resisting  moment  -  in-lb 

ON 

Normal  force  at  rotating  band,  NQ  -  lb 

ONM,  ONN 

Maximum,  minimum  NQ  -  lb 

ONMT,  ONNT 

Time  when  NQ  is  maximum,  minimum  -  ms 

ONMX,  ONNX 

Projectile  travel  when  NQ  is  maximira,  minimum  -  in 

ONRM 

Nq  expressed  in  metric  unit  -  N 

OPSLP 

Slope  of  gun  tube-rotating  band  contact  pressure  data 
curve 

P 

Total  base  pressure  —  lb 

PAD,  PFD 

Distance  of  axial  point  A,  F  to  c.g.  -  in 

PCGX 

Distance  of  base  pressure  center  to  c.g.  -  in 

PM,  PN 

Maximum,  minimum  P  -  lb 

PMT,  PNT 

Time  when  P  is  maximum,  minimum  -  ms 

PMX ,  PNX 

Projectile  travel  when  P  is  maximim,  minimum  -  in 

PRM 

Moment  of  total  base  pressure  about  c.g.  -  in-lb 

PROJRBD 

Diameter  of  rotating  band  seat  -  in 

PWT 

Weight  of  projectile  -  lb 

PX,  PY 

X,  Y  component  of  P  -  lb 

R,  RO 

Radius  of  gun  tube  bore  -  in 

RB 

Radius  of  bourrelet  -  in 
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Table  2.  (cont) 


Symbol  Description 


RBANG ,  RBM 

Turning  angle,  deg;  moment,  in-lb  of  rotating  band  trans¬ 
verse  resisting  moment  data 

RBD 

Outer  radius  of  rotating  band  when  projectile  is  loaded 
in  gun  tube  -  in 

RBDEF,  RBLOAD 

Deflection,  in;  applied  load,  lb  of  rotating  band  com¬ 
pression  test  data 

RBE 

Young’s  modulus  of  rotating  band  -  10^  psi 

RBID,  RBOD 

Inner,  outer  diameter  of  rotating  bnad  before  it  is  fit¬ 
ted  on  projectile  band  seat  -  in 

RBK,  RBMK 

Spring  constant  -10  lb/in;  angular  spring  constant  of 
rotating  band  -  10°  in-lb/deg 

RBKM 

RBK  expressed  in  metric  unit  —  10^  N/cm 

RBNU 

Poisson's  ratio  of  rotating  band  material 

RBPMX,  RBPR 

Distance,  in;  pressure,  psi  of  gun  tube-rotating  band 
pressure  input  data 

RBPRAV 

Average  gun  tube-rotating  band  contact  pressure  -  psi 

RBPRC 

Gun  tube-rotating  band  contact  pressure  at  breech  end, 
computed  in  the  program  -  psi 

RBPRCM 

RBPRC  expressed  in  metric  units  -  kPa 

RBPRMR 

Ratio  of  gun  tube-rotating  band  contact  pressure  at  muz¬ 
zle  and  to  that  at  breech 

RBPW 

Effective  width  of  rotating  band  used  to  compute  total 
gun  tube-rotating  band  contact  pressure  -  in 

RBW 

Nominal  width  of  rotating  band  -  in 

RBX 

Distance  from  projectile  base  to  nearest  side  of  rotating 
band  -  in 

RHO 

Air  density  used  to  compute  aerodynamic  forces  (0.002378 
slugs/ft ^  for  default  value) 

RTBAND 

Control  number  to  identify  rotating  band  or  obturator 
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Table  2.  (cont) 


Symbol  Description 


STHETA 

Sine  of  angle  alpha 

TADDMX,  TADDNX 

Projectile  travel  when  yaw  acceleration  is  maximum,  mini¬ 
mum  -  in 

TADMX,  TAD NX 

Projectile  travel  when  yaw  velocity  is  maximum,  minimum  - 
in 

TAMX,  TANX 

Projectile  travel  when  yaw  angle  is  maximum,  minimum  -  in 

TBRB,  TBRO 

Radius  of  worn  gun  tube  bore  at  location  of  bourrelet 
rotating  band  -  in  9 

THETA 

Yaw  angle  -  rad 

THETAC 

Yaw  angle  when  bourrelet  just  contacts  gun  tube  at  6 
o'clock  position  -  rad 

THETACO 

THETAC  expressed  in  degree 

THETAD 

Yaw  velocity  -  rad/s 

THETAM,  THETAN 

Maximum,  minimum  yaw  angle  —  rad 

THETAMT,  THETANT 

Time  when  yaw  angle  is  maximum,  minimum  -  ms 

THETAO 

THETA  expressed  in  degree 

THETAPM 

THETA  expressed  in  milliradians 

THE TAR 

Control  number  in  using  computed  contact  yaw  angle  as 
initial  position  8 

THTADD 

Yaw  acceleration  -  rad/s2 

THTADM,  THTADN 

Maximum,  minimum  yaw  velocity  -  rad/s 

THTADMT,  THTADNT 

Time  when  yaw  velocity  is  maximum,  minimum  -  ms 

THTASGN 

Sign  of  THETA  or  yaw  angle 

TI 

Time  step  size  used  in  computation  -  sec 

TIAV 

Average  time  step  size  -  sec 

TIJ 

Cncreraent  of  time  step  size  used  in  variable  time  step 
size  computation  -  sec 
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Table  2.  (cont) 


Symbol  Description 


TIME 

Time  in  computation  -  ms 

TITLE 

Subtitle  of  output  figure 

TITLEM 

Title  of  analysis  or  computer  run,  also  used  as  figure 
subtitle  if  no  subtitle  is  given 

TITLES 

Subtitle  of  analysis  or  computer  run,  also  used  as  figure 
subtitle  if  no  subtitle  is  given 

TMAX 

Estimated  or  desired  time  of  computation,  also  used  as 
max  value  of  time  scale  in  plotting  figures  -  ms 

TMOI 

Transverse  moment  of  inertia  of  projectile  -  lb-in2 

TRAVEL 

Travel  or  maximum  displacement  of  projectile  -  in 

TRAVELM 

TRAVEL  expressed  in  metric  unit  -  m 

TSLAST 

Last  time  step  size  of  uniformly  decreased  step  size 
computation  -  ms 

TSTEP 

Time  step  size  for  constant  step  size  computation  or 
first  step  size  for  uniformly  decreased  step  size  compu¬ 
tation  -  ms 

TTADDM,  TTADDN 

Maximum,  minimum  yaw  acceleration  -  rad/s^ 

TTADDMT,  TTADDNT 

Time  when  yaw  acceleration  is  maximum,  minimum  -  ms 

TTHETA 

Tangent  of  angle  theta 

TWDC,  TWMX 

Change  in  gun  tube  bore  diameter,  mil;  distance,  in  of 
worn  tube  input  data 

VB 

Velocity  of  bourrelet  contact  point  where  normal  force 
resultant  acts  -  in/s 

VBSGN 

Sign  of  velocity  VB 

VO 

Velocity  of  rotating  band  contact  point  where  normal 
force  resultant  acts  -  in/s 

VOSGN 

Sign  of  velocity  VO 

WRB,  WRO 

Change  in  radius  of  worn  gun  tube  at  location  of  bourre¬ 
let,  rotating  band  -  in 
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Table  2,  (cont) 


-Symbol  Description 


WROSLP 

Slope  of  worn  gun  tube  data  curve 

X,  XB 

Travel  or  axial  displacement  of  projectile  c.g.,  bourre- 
let  -  in 

XBMZL 

Bourrelet  distance  to  muzzle  —  in 

XD,  XDD 

Axial  projectile  velocity  -  in/s;  acceleration  -  in/s2 

XDDM,  XDDN 

Maximum,  minimum  XDD  -  in/s2 

XDDMT,  XDDNT 

Time  when  XDD  is  maximum,  minimum  -  ms 

XDDMX,  XDDNX 

Projectile  travel  when  XDD  is  maximum,  minimum  —  in 

XDDR 

XDD  expressed  in  g's 

XDR 

XD  expressed  in  ft/s 

XDRM 

XC  expressed  in  m/s 

XI,  YI 

X,  Y  value  of  initial  c.g.  location  -  in 

XIM,  YIM 

XI,  YI  expressed  in  metric  unit  -  cm 

XL 

Distance  of  c.g.  to  mid-section  of  rotating  band  -  in 

XMASS 

Mass  of  projectile  -  slug/12 

XRM 

Projectile  axial  displacement  or  travel  expressed  in 
metric  unit  -  cm 

XTMOI 

TM0I  expressed  in  slug-in2/12 

Y 

Normal  displacement  of  projectile  c.g.  -  in 

YD,  YDD 

Normal  velocity  -  in/s,  acceleration  of  projectile  c.g.  - 
in/sz  6 

YDDM,  YDDN 

Maximum,  minimum  YDD  -  in/s2 

YDDMT ,  YDDNT 

Time  when  YDD  is  maximum,  minimum  -  ms 

YDDMX,  YDDNX 

Projectile  travel  when  YDD  is  maximum,  minimum  -  in 

YDDR 

YDD  expressed  in  g's 
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Table  2*  (cont) 


Symbol  Description 


YDM,  YDN 

Maximum,  minimum  YD  -  in/s 

YDMT,  YD NT 

Time  when  YD  is  maximum,  minimum  -  ms 

YDMX,  YDNX 

Projectile  travel  when  YD  is  maximum,  minimum  -  in 

YDR 

YD  expressed  in  ft/s 

YDRM 

YD  expressed  in  m/s 

YM,  YN 

Maximum,  minimum  Y  -  in 

YMT,  YNT 

Time  when  Y  is  maximum,  minimum  -  ms 

YMX,  YNX 

Projectile  travel  when  Y  is  maximum,  minimum  -  in 

YRM 

Y  expressed  in  metric  unit  -  cm 

ZM 

Total  moment  about  Z  axis  -  in-lb 
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Table  3.  Definitions  of  symbols  used  in  subroutine  CLIPV 


Symbol 


XD,  YD 

XD(N),  YD(N) 

XO 

X 

YO 

SLP 

ND 


Description 
Variable  or  data  pairs 
Nth  variables  or  data  pairs 
Known  intermediate  value  of  XD 
Absolute  value  of  XO 

Intermediate  value  of  YD  to  be  computed 
Slope  of  XD-YD  curve  in  computation 
Number  of  data  pairs 
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Table  4.  Definitions  of  symbols  used  in  subroutine  CMAXMIN 


Symbol  Description 


X 

Variable,  its  maximum  and  minimum  values  are  to  be  deter¬ 
mined 

XM 

Maximum  value  of  X 

XN 

Minimum  value  of  X 

TM,  TN 

Time  when  maximum,  minimum  value  of  X  is  first  computed 

DM,  ON 

Axial  projectile  displacement  when  maximum,  minimum  value 
of  X  is  first  computed 

T 

Time  corresponding  to  X 

D 

Axial  projectile  displacement  corresponding  to  X 
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Table  5.  Definitions  of  symbols  used  in  subroutine  DATAINP 


Symbol 

Description 

ACM 

Axial  distance  of  projectile  c.g.  to  aerodynamic  center, 
expressed  in  metric  unit  -  cm 

ALPHAO 

Alpha  expressed  in  degrees 

APAX,  APFX 

Distance  of  projectile  axial  point  A,  F  to  projectile 
nose  -  in 

APAXM 

APAX  expressed  in  metric  unit  -  cm 

BET  AO 

BETA  expressed  in  degrees 

BRTEM 

BRTE  expressed  in  metric  unit  -  10^  kPa 

BRTID 

Inside  diameter  at  bourrelet  -  in 

BRTIDM 

BRTID  expressed  in  metric  unit  -  cm 

BRTTKM 

BRTTK  expressed  in  metric  unit  -  era 

BRTWM 

BRTW  expressed  in  metric  unit  -  cm 

BRTXM 

BRTX  expressed  in  metric  unit  -  cm 

CGD 

Projectile  diameter  at  c.g.  -  in 

CGDM 

CGD  expressed  in  metric  unit  -  cm 

CHORDM 

CHORD  expressed  in  metric  unit  -  cm 

CWTM 

CWT  expressed  in  metric  unit  -  kg 

FAFMN,  FAFMX 

Minimum,  maximum  figure  scale  of  aerodynamic  drag  force  - 
kN 

FAPFMN,  FAPFMX 

Minimum,  maximum  figure  scale  of  normal  acceleration  at 
axial  point  F  -  gfs 

FCGAAMN,  FCGAAMX 

Minimum,  maximum  figure  scale  of  axial  projectile  accel¬ 
eration  -  gTs 

FCGADMN,  FCGADMX 

Minimum,  maximum  figure  scale  of  axial  projectile  dis¬ 
placement  -  m 

FCGAVMN,  FCGAVMX 

Minimum,  maximum  figure  scale  of  axial  projectile  velo¬ 
city,  m/s 
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Table  5*  (cont) 


Symbol  Description 


FCGNAMN,  FCGNAMX 

Minimum,  maximum  figue  scale  of  normal  projectile  accel¬ 
eration,  gTs 

FCGNDMN,  FCGNDMX 

Minimum,  maximum  figure  scale  of  normal  projectile  dis¬ 
placement  -  mm 

FCGNVMN,  FCGNVMX 

Minimum,  maximum  figure  scale  of  normal  projectile  velo¬ 
city  -  m/s 

FGTMAX,  FGXMAX 

Maximum  time,  ms;  travel  (or  axial  displacement)  scale  of 
output  figure  -  in 

FGXMAXM 

FGXMAX  expressed  in  metric  unit  -  m 

FLFMN,  FLFMX 

Minimum,  maximum  figure  scale  of  lateral  or  normal  force 
-  kN 

FPEAM 

0 

FPEA  expressed  in  metric  unit  -  cm 

FTMAX,  FTMIN 

Maximum,  minimum  figure  scale  of  time  -  ms 

FXMAX,  FXMIN 

Maximum,  minimum  figure  scale  of  travel  or  axial  dis¬ 
placement  -  in 

FYAWAMN ,  FYAWAMX 

Minimum,  maximum  figure  scale  of  yaw  acceleration  - 
rad/s ^ 

FYAWDMN,  FYAWDMX 

Minimum,  maximum  figure  scale  of  yaw  angle  -  milliradians 

FYAWVMN,  FYAWVMX 

Minimum,  maximum  figure  scale  of  yaw  velocity  -  rad/s 

MXMN 

Control  number  of  figure  scale  input  data 

NAF,  NAP FA,  NCGDVAA, 
NCGDVAN,  NLF ,  NYAWA, 
NYAWDV 

Control  number  of  figure  plotting  of  air  resistance 
force,  normal  acceleration  at  axial  point  F,  axial 
projectile  motion,  normal  projectile  motion,  lateral 
forces,  yaw  acceleration,  yaw  angle,  and  velocity 

NFGVSX 

Control  number  for  plotting  curves  against  travel 

NPLOT 

General  control  number  of  figure  plotting 

PCGXM 

PCGX  expressed  in  metric  unit  -  cm 

PROJL 

Length  of  projectile  -  in 

PROJLM 

PROJL  expressed  in  metric  unit  -  cm 
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Table  5*  (cont) 


Symbol 

PWTM 

RBDM 

RBEM 

RBIDM 

RBMKM 

RBODM 

RBPRAVM 


Description 

PWT  expressed  in  metric  unit  -  kg 

RBD  expressed  in  metric  unit  -  cm 

RBE  expressed  in  metric  unit  -  10^  kPa 

RBID  expressed  in  metric  unit  -  cm 

RBMK  expressed  in  metric  unit  -  cm-N/rad 

RBOD  expressed  in  metric  unit  -  cm 

RBPRAV  expressed  in  metric  unit  -  kPa 


RBPWM,  RBWM,  RBXM  RBPW,  RBW,  RBX  epxressed  in  metric  unit  -  cm 

RflOM  RHO  expressed  in  metric  unit  -  g/cm3 

IJNITAF,  IJNITBK,  UNITFP  Control  number  of  input  data  units  for  input  data  of 
IJNITMS,  UNITRBK,  aerodynamic  forces,  bourrelet  elasticity,  firing  pres- 

IJNITRBP ,  UNITTW  sure,  dimensions  and  physical  properties,  band  elasti¬ 

city,  gun  tube-band  contact  pressure,  gun  tube  wear 
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Table  6.  Definitions  of  symbols  used  in  subroutine  DTARO 


Symbol  Description 


X,Y 

Data  pairs 

NODATA 

Number  of  data  pairs 

XO 

New  reference  point  or  origin  location  (distance  from 
first  data  pair  to  new  origin) 

YRTO 

Ratio  of  converted  Y  value  to  original  Y  value 
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Table  7. 

Definitions  of  symbols  used  in  subroutine  DTATBL 

Symbol 

Description 

X(I),  Y(I) 

Given  data  pairs,  X,Y 

XA,  YA;  XB,  YB 

Converted  data  pairs  to  be  printed  in  one  line 

M,  NFORMAT 

Format  statement  reference  number 

NODATA 

Number  of  data  pairs 

XRTO ,  YRTO 

Conversion  ratio  of  X,Y 
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Table  8.  Definitions  of  symbols  used  in  subroutine  FIGCTRL 


Symbol 

APFSTL 

ENDLBL 

FDDMM,  FDDNM 
FDMM,  FDNM 
FNMM,  FNNM 
FTS 

GEE 

HBEDGE 

HEL ,  HFT ,  HL ,  HTL 

MXN 

NFG 

NOD A,  NODTA 
PFDM 

TD2,  TD3,  TD4,  TD5, 
TD6,  TD7,  TD8,  TD9, 
TD10,  TD11,  TD1 5, 
TD16,  TD18,  TD31 

THETAMM,  THETANM 

TLABEL 

TS 

WHR 

XDDMM,  XDDNM 


Description 


Figure  subtitle  indicating  location  of  axial  point  F 
End  label  of  output  curve 
FDDM,  FDDN  expressed  in  g's 

Maximum,  minimum  aerodynamic  drag  force  expressed  in  kN 

Maximum,  minimum  lateral  force  expressed  in  kN 

Space  between  first  figure  subtitle  line  and  top  of  fig¬ 
ure  -  in 

Acceleration  unit  conversion  factor 

Distance  of  bottom  edge  of  graph  paper  to  X  axis  of  fig¬ 
ure  -  in 

Height  of  character  of  curve  end  label,  figure  title, 
axial  point  F  distance  subtitle,  figure  subtitle  -  in 

Control  number  =  MXMN 

Graph  location  control  number 

Number  of  curve  data 

PFD  expressed  in  metric  unit  -  cm 

Tape  recorded  data  for  time,  travel,  normal  projectile 
displacement,  yaw  angle,  axial  projectile  velocity, 
normal  projectile  velocity,  yaw  velocity,  axial  pro¬ 
jectile  acceleration,  normal  projectile  acceleration,  yaw 
acceleration,  lateral  force  at  rotating  band,  lateral 
force  at  bourrelet,  aerodynamic  drag  force,  normal  accel¬ 
eration  at  axial  point  F 

THETAM,  THETAN  expressed  in  milliradians 
Time  axis  label 

Space  between  figure  subtitle  lines  -  in 
Ratio  of  label  character  width  to  height 
XDDM,  XDDN  expressed  in  g’s 
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Table  8.  (cont) 

Symbol 

Description 

XDMM,  XDNM 

XDM,  XDN  expressed  in  metric  unit  -  m/s 

XLABEL,  YLABEL 

Label  of  figure  axis  X,Y 

XMM,  XNM 

Maximum,  minimum  X  expressed  in  metric  unit  -  m 

XSIZE,  YSIZE 

Length  of  figure  axis  X,Y  -  in 

XSPACE,  YSPACE 

Distance  between  figurs  in  X,Y  direction  -  in 

YDDMM,  YDDNM 

YDDM,  YDDN  expressed  in  gfs 

YDMM,  YDNM 

YDM,  YDN  expressed  in  metric  unit  -  m/s 

YMM,  YNM 

YM,  YN  expressed  in  metric  unit  -  mm 
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Table  9.  Definitions  of  symbols  used  in  subroutine  FGCV3 


Symbol 
ELX,  ELY 
FGTITLE 
NOCRV 
NPLNO 
OX,  OY 
X 

XMAX,  XMIN; 
YMAX,  YMIN 

XSTART,  YSTART 

XT,  YT 

YA,  YB,  YC 


Description 

Starting  X,Y  coordinates  to  print  curve  end  label 
Title  of  figure 

Number  of  curves  plotted  in  one  graph 
Number  of  figure  plotting 
X,Y  coordinates  of  figure  origin 
Common  X  coordinate  data 
Maximum,  minimum  of  X;Y  scale 

Starting  value  of  X,Y  scale  plotting 

Starting  X,Y  coordinates  to  write  figure  title  or  sub¬ 
title 

First,  second,  third  curve  Y  coordinate  data 
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Table  10. 

Definitions  of  symbols  used  in  subroutine  FGDA3 

Symbol 

Description 

TDl ,  TD2,  ...  TD31 

Data  recorded  on  tape  1 

X(I),  Y(I) ;  X(I), 
YB(I);  X(l),  YC( I) 

Data  to  be  used  to  plot  curve  1,  2,  3 

XD,  YDA,  YDB,  YDC 

X( I) ,  YA( I) ,  YB( I) ,  YC( I)  data  read  from  tape  1 

XDR,  YDAR,  YDBR,  YDCR 

Conversion  factor  of  XD,  YDA,  YDB,  YDC 
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Table  11.  Definitions  of  symbols  used  in  subroutine  FGMN3 


Symbol  Description 


ABSDA,  ABSDB,  ABSDC 

Absolute  value  of  difference  of  maximum  and  minimum  val¬ 
ues  of  Y  coordinates  for  curve  1,  2,  3 

CXM,  CXN 

Maximum,  minimum  computed  value  of  X 

CYAM,  CYAN;  CYBM, 

CYBN;  CYCM,  CYCN 

Maximum,  minimum  computed  value  of  Y  for  curve  1,  2,  3 

NGNB ,  NGNC 

Control  number  used  in  computing  scale  ratio  for  2,  3 

curves 

NYDRBA,  NYDRCA 

Integer  ratio  of  ABSDB  to  ABSDA,  ABSDC  to  ABSDA 

YAM,  YAN;  YBM,  YBN; 
YCM,  YCN 

Maximum,  minimum  YA,  YB,  YC 

YBMN,  YBMX;  YCMN, 

Converted  minimum,  maximum  YB;  YC  used  in  scale 

YCMX  computation 

YDRBA,  YDRCA 

Ratio  of  ABSDB  to  ABSDA,  ABSDC  to  ABSDA 

YYLABEL 

Y  axis  label  for  2  or  3  curves 
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Table  12. 

Definitions  of  symbols  used  in  subroutine  FGPRTCV 

Symbol 

Description 

CURVE 

Array  name  of  curve  symbols  used  to  print  curves 

DX,  DY 

Number  of  units  per  inch  of  X,  Y  axis 

DXIO,  DYIO 

Number  of  units  per  line  in  X,  Y  direction 

FIELD 

Field  descriptor  for  printer  plotting  scale  printout 

FMTO,  FMT1,  FMT2, 
FMT3,  FMT4,  FMT5 

Format  number  for  printer  plotting  label  formats 

N 

Number  of  data 

NTX,  NTXX 

Number  of  spaces  in  printer  plotting  title  and  subtitle 
printout 

NX,  NY 

X,  Y  number  of  symbols  used  in  printer  plotting 

NXX 

Number  of  X  scale  values 

XSN 

X  scale  values 
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Table 

13.  Definitions  of  symbols  used  in  subroutine  FGPRTSB 

Symbol 

Description 

BAR 

Symbol  of  curve,  — 

BLANK 

Symbol  of  curve  for  no  plotting 

CVSBL 

Array  name  of  symbols  of  different  printer  curves; 
CVSBL(l)  =  A,  CVSBL( 2)  =  B,  CVSBL(3)  =  C 

L 

Line  number 

LN 

Scale-annotation  line  number 

LX 

Symbol  position  index  number 

NC 

Curve  index  number 

STAR 

Curve  symbol  determined  from  CVSBL 

VLINE 

Curve  symbol  for  vertical  line 
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Table  14,  Definitions  of  symbols  used  in  subroutine  FGSTRT 


Symbol  Description 


D,  DX 

Units  of  scale  per  unit  length  of  scale 

MX 

Flag  number  to  indicate  whether  maximum  and  minimum  val¬ 
ues  are  given  (when  MX  =  1)  or  computed  values  in  the 
program  (MX  =  0)  are  used 

XMAX,  XX( 2) 

Maximum  value  of  scale 

XMIN,  XX(1) 

Minimum  value  of  scale 

XSIZE,  SIZE 

Size  or  length  of  axis 

XSTART,  START 

Starting  value  of  axis  scale 
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Table  15. 

Definitions  of  symbols  used  in  subroutine  OPDTA 

Symbol 

Description 

DA,  DB,  DC,  DD,  DE, 
DF,  DO,  DH,  DI,  DJ 

Intermediate  symbols  for  coverting  English  units  to 
metric  units  of  data 

NO 

Data  number 
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Figure  4  (Continued) 
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Figure  4  (Continued) 
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Flowchart  of  Subroutine  DIPTBL 
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Flowchart  of  Subroutine  FGCV3 
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APPENDIX  A 
INPUT  DATA  DECK 
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The  preparation  of  the  input  data  cards  has  been  described  in  section: 
Input  Data  Coding.  To  facilitate  preparation  of  input  data  cards,  and  for  easy 
reference,  the  input  data  deck  is  tabulated  below. 


Input  data 

deck 

Card 

(set) 

order 

Data  type 

Format 

List  of  sybmols 

i 

Title 

5A10 

TITLEM(I),  1=1,5 

2 

Subtitle 

5A10 

TITLES(I),  1=1,5 

3 

Computation, 
output  control 

7110 

NDAMS ,  NDARAD , NIDTBL ,  NCCK, 

NOPTBL , NOOPDTA ,  NPOPDTA 

4 

Time 

3F10.5 

TMAX,  TSTEP,  TSLAST 

5 

Projectile 

(A) 

8F10.5 

PRO JL ,  RBX,  RBW,  RBD,  CGX,  CGD ,  CGEL 
CGEMT 

(B) 

7F10.5 

BRTX,  BRTW,  BRTOD ,  BRTID ,  BRTTK , 

APAX,  APFX 

(C) 

2F10.4 

PWT,  TMOI 

6 

Gun  tube 

4F10.5 

ALPHA,  BORED,  FPEA,  TRAVEL 

7 

Projectile 

inclination 

3F10.5 

BETA,  THETA,  THETAR 

8 

Figure 

(A) 

8110 

NCGDVAA,  NCGDVAN,  NLF ,  NAP FA,  NAF, 
NYAWDV ,  NYAWA,  NFGVSX 

Omit  following  (B)  -  (N) 

cards  if  no  figures  desired. 

(B) 

110,  3F10.5 

MXMN,  FGXMAXM,  FGTMAX,  FGXMAX 

Omit 

following  (C) 

-  (G) 

cards  if  MXMN  =  0. 

(C) 

5F10.4 

FTMAX,  FTMIN,  FXMAX,  FXMIN,  UNITMS 

(D) 

6F10.4 

FCGADMX,  FCGADMN,  FCGAVMX ,  FCGAVMN, 
FCGAAMX,  FCGAAMN 

(E) 

6F10.4 

FCGNDMX,  FCGNDMN,  FCGNVMX,  FCGNVMN, 
FCGNAMX,  FCGNAMN 

(F) 

6F10.4 

FLFMX,  FLFMN ,  FAPFMX,  FAPFMN,  FAFMX, 
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(G) 

6F10.4 

FYAWDMX ,  FYAWDMN ,  FYAWVMX ,  FYAWVMN 
FYAWAMX ,  FYAWAMN 

(H) 

6110 

NT( I) ,  1=1,6 

(I) 

5A10 

TITLE( J, 1 ) ,  J=1 ,5 

If  NT(1)  =  0,  omit  this  card. 

(J) 

5A10 

TITLE( J,2) ,  J=1 ,5 

If  NT(2)  =  0,  omit  this  card. 

(K) 

5A10 

TITLE( J,3) ,  J=1 ,5 

If  NT(3)  =  0,  omit  this  card. 

(L) 

5A10 

TITLE(J,4) ,  J-1,5 

If  NT(4)  =  0,  omit  this  card. 

(M) 

5A10 

TITLE(J,5) ,  J=1 ,5 

If  NT(5)  =  0,  omit  this  card. 

00 

3A10 

TITLE( J ,6 ) ,  J=1 ,5 

If  NT(6)  =0,  omit  this  card. 

9 

Firing  pressure 
(A) 

2110, 

4F10.5 

NOFP,  NCHAMP,  CWT,  PCGX,  UNITFP 

(B) 

4(F10 .5 , 
F10.2) 

FPT(I),  FP( I) ,  1=1,  NOFP 

10 

Tube  wear  (A) 

110,  F10.5 

NOTW,  UNITTW 

(B) 

8F10.5 

TWMX(I),  TWDC(I),  1=1,  NOTW 

If  NOTW  =  0,  omit  this  card(s). 

11 

Bourrelet 
elasticity  (A) 

110,  3F10.4 

NOEBK,  BRTK,  BRTE,  UNITBK 

(B) 

4(F10 .5 , 
F10.2) 

BDEF(I),  BLOAD(I),  1=1,  NOEBK 

If  NOEBK  =  0,  omit  this  card(s). 

12 

Band  elasticity 
(A) 

2110, 

4F10.4 

NOERBK,  NOERBMK,  RBK,  RBMK,  RTBAND 
UNITRBK 

(B) 

4(F10 .5 , 
F10.2) 

RBDEF(I),  RBLOAD(I),  1=1,  NOERBK 

If  NOERBK  =  0,  omit  this  card(s). 

(c) 

4(F10 .5 , 
F10.2) 

RBANG(I),  RBM(I),  1=1,  NOERBMK 

If  NOERBMK  =  0,  omit  this  card  (s) 

13 

Friction 

(A) 

110,  F10.2, 
6F10.5 

NORBPR,  RBPRAV,  RBMU ,  BRTMU ,  RBPW, 
RBPRMR,  RBNU,  UNITRBP 
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(B) 

4F10.5 

RBOD,  RBID,  RBE,  PROJRBD 

If  RBPRAV  is  not  zero  or  RBNU  =  0, 
omit  this  card. 

(C) 

4(F10.5, 

F10.2) 

RBPMX(I),  RBPR(I),  1=1,  NORBPR 

If  RBPRAV  is  not  zero  or  NORBPR  =  0 
omit  this  card(s). 

14 

Air  resistance  (A) 

110,  5F10.5 

NOCD,  CDAV,  RHO,  AC,  DTALAM,  UNITAF 

(B) 

8F10.5 

CDANG(I) ,  CD( I) ,  1=1,  NOCD 

If  NOCD  =  0,  omit  this  card(s). 

(C) 

2110, 

3F10.5 

NOCL,  NOCM,  CLAV,  CMAV,  CHORD 

If  DTALAM  =  0,  omit  this  and  the 
following  (D)  and  (E)  cards. 

(D) 

8F10.5 

CLANG(I) ,  CL( I) ,  1=1,  NOCL 

If  NOCL  =  0,  omit  this  card(s). 

(E) 

8F10.5 

CMANG(I) ,  CM( I) ,  1=1,  NOCM 

If  NOCM  =  0,  omit  this  card(s). 

15 

End  of  data 

7/8/9  in  column  1 
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APPENDIX  B 

COMPUTER  OUTPUT  OF  EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALLOT2 


Printer  output,  pages  123  -  149 
Plotter  output,  pages  150  -  163 
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BALLOT 2  IN BORE PROJECT  I LEMQT I OMAN A  L Y£ I SPRGGRAMBYSHCHU ASDLCWS  LAP  RADCGMDOVERN JV  ERS I ONONE 
BBBB33SBBDB jAAAAAAAAAAAA L L LLLLL LL ILL L L L L LlL L L L LLOGOGOCQOGOOOT I T T T IT T T7TT2 222222 22222 
BB3BBBBB&B33AAAAAAAAAAAALLLLLLLLLLLLLLLLLLL  LLLLlOOOuOOOOOOOOTTT  TTTTTTTTT2  22222222222 


3B 

B2BAAIAAA  AAAAALL 

LUUll'.LL 

LLLLLLLLOOO 

OOCTT 

TT222 

222 

B8 

BBBB 

bba^aa 

AAAALL 

I.LLLLlI-lLL 

LLLLLLLLOO 

0000 

OOTTTTT 

TTTTT22 

2222  22 

B3 

B2BB 

BBAAA 

AA 

AAALL 

LLLLLLLLLL 

LLLLLLlLOO 

0000 

OOTTTTT 

TTTTT22222222  22 

BB 

3BBAA 

AAAA 

AALL 

LLLLL^liLL 

lLLLLLLLOO 

0000 

OOTTTTT 

TTTTT2  222 

22 

BB 

BBBB 

bbaa 

aaaa 

AALL 

LLLLLLLLLL 

LLLLLLLLC'O 

0000 

OOTTTTT 

TTTTT222 

222 

BB 

BB  BB 

bbaa 

AALL 

Ul ILL  LlLL 

LLLLLLLLOO 

cooo 

OOTTTTT 

TTTTT22 

2222222 

B3 

BBBB 

BB  A  A 

AAAA 

AALL 

LLLLLLLLLL 

LLLLLLLLOO 

0000 

OOTTTTT 

TTTTT22 

22222222 

BB 

BBBB 

03  A  A 

aaaa 

AALL 

LLLLLLLLLL 

LLLLLLLLOO 

0000 

OOTTTTT 

TTTTT22 

22222222 

BB 

38B  A  A 

AAAA 

AALL 

LLLL 

L  1.000 

OOQTTTTT 

TTT7T222 

22 

BS33SBB3B3BBAAAAAaAAAAAAlLLLLLLLL  L  LLLLLLLLLL LLLLOOOOOOOOOOOO" TTTTTTTTTTT22222 22 22222 
BBBBBBBBBBBBAAAAAaAAAAAALLLLlLLLLLLLLLLLLLLLLLLLOOOOCOOOOOOOTTT TTTTTTTTT222222222222 
BALLOT 2  IN BORE PROJECT  I LEMOT IGNANALYS I SPROGRAMBYSHCHU ASDLCWS LARRADCOMDOVERN JVERSIONONE 


IN-SORE  MOTION  ANALYSIS 
CF  PROJECTILE  IN  STATIONARY  SMOOTH  GUN  TUBE 


X AMPLE  OF  COMMUTATION  OF  PROGRAM  BA  L  LOT  2 


fable 


1 


***** IN-DORE  MOTION  ANALYSIS  input  data  cards***** 

t 


12345678901  23-1567890  12.34  567890  12345678901234567890 1234  567890  123456  7890 12  34567890 


(TITLE) 

EXAMPLE  OF  .COMPUTATION  OF  PROGRAM  BALLOT2 
(SUBTITLE) 


(COMPUTATION  AND  OUTPUT  CONTROLS) 


0 

0 

0 

0 

0 

0 

20 

(TIME  DATA) 

6.00000 

.00050 

0.00000 

( PROJECTILE 

DATA) 

23.04900 

8.56500 

.81000 

0.00000 

1  1  .280^0 

3.17600 

-.01000 

0.00000 

13.60500 

.62500 

4.71500 

0.00000 

0.00000 

5.00000 

2.50000 

15.5977 

7.2095 

( GUN  TUBE  DATA) 

0.00000 

4 . 73400 

0.00000 

187.00000 

(PROJECTILE 

.  INCLINATION) 

0.00000 

0 .0000000 

1  .00000 

(OUTPUT  FIGURE  DATA) 

1 

1 

1 

1 

1 

1 

1 

1 

G 

4 .80000 

6.00000 

0.00000 

(OUTPUT  FIGURE  SUBTITLE  CONTROLS  AND  SUBTITLES) 

G 

0 

0 

0 

0 

0 

(TIME-PRESSURE  DATA) 

63 

0 

.25000 

0.00000 

0.00000 

0.00000 

100.00 

.  10000 

264.00 

.20000 

443.00 

.30000 

676.00 

.40000 

979.00 

.50000 

1371 .00 

.60000 

1877.00 

.70000 

2515.00 

.80000 

3335.00 

.90000 

4370.00 

1  .  Q0000 

5670.00 

1 .1  0000 

7288.00 

1 . 20000 

9283.00 

1 .30000 

1 1669.00 

1  .41600 

14746.00 

1  .41605 

14697.00 

1 .50000 

1 7267.00 

1 .60000 

20754.00 

1  .70000 

24770.00 

1  .80000 

29379.00 

1 .90000 

34293.00 

2.00000 

39285.00 

2. 1O0CO 

44302.00 

2.20000 

48974.00 

2.30000 

53053.00 

2.40000 

56337.00 

2.50000 

58704.00 

2.60000 

601 18.00 

2.70000 

60630.00 

2.80000 

60348.00 

2.90000 

59416.00 

3.00000 

57973.00 

3. 1 0000 

56209.00 

3.20000 

54160.00 

3.30000 

51 975.00 

3.40000 

49731 . 00 

3'.  50000 

47487.00 

3.60000 

45284.00 

3.7000C 

42650.00 

3.80000 

39651.00 

3.90000 

36687.00 

4 . 0000C 

33869.00 

4. 10000 

31 248.00 

4.20000 

28839.00 

4.30000 

26656.00 

4 . 40000 

24688.00 

4.50000 

22913.00 

4.60000 

21310.00 

4.70000 

19861 .00 

4.B0000 

18549.00 

4.900C0 

17357.00 

5.00000 

16274.00 

5. 1 0000 

15287.00 

5. 10500 

15322.00 

5.20000 

14402.00 

5.30000 

13565.00 

5.40000 

12802.00 

5.50000 

12104.00 

5.60000 

11463.00 

5.70000 

10874.00 

s.eoooo 

10331 .00 

5.90000 

9828.00 

5.94500 

9614.00 

(TUBE  WEAR 

DATA) 

0 

0.00000 

(BOURRELET 

ELASTICITY 

DATA) 

0 

. 28000 

0.00000 

0.00000 

(BAND  ELASTICITY  DATA) 

0 

0 

.83000 

0.00000 

0.00000 

0.00000 

(FRICTION  DATA) 

0 

5125.00 

.04000 

.1 5000 

0.00000 

0.00000 

0.00000 

0.00000 

(AERODYNAMIC  COEFFICIENT  DATA) 

- 

0 

1 .00000 

0.00000 

0.00000 

0.00000 

0.00000 

( END) 
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TABLE  2  PROJECTILE  AND  GUN  TUBE  DATA 


EXAMPLE  OF  COMPUTATION  OF 

PROGRAM  BALLOT2 

DIMENSIONS  AND  PROPERTIES 

OF  PROJECTILE 

LENGTH  OF  PROJECTILE 

(  23. 049  IN) 

58.544  CM 

OBTURATOR  TO  BASE  DISTANCE 

OBTURATOR  WIDTH 

(  8.565  IN) 

(  .BIO  IN) 

21.755  CM 
2.057  CM 

C.  G.  TO  BASE  DISTANCE 

C.  G.  ECCENTRICITY  (DISTANCE) 

C.  G.  ECCENTRICITY  (MOMENT) 

•(  11.201  IN) 

(  -.0100  IN) 
(-2.5  IN-OZ) 

28.652  CM 
-.254  MM 
-17.62  MM-N 

BOURRELET  TO  BASE  DISTANCE 

BOURRELET  DIAMETER 

(  13.605  IN) 

(  4.71 5  IN) 

34.557  CM 

1 1.976  CM 

BOURRELET  SPRING  CONSTANT  (10E6) 
OBTURATOR  SPRING  CONSTANT  (10E6) 
AVERAGE  OBTURATOR-TUBE  PRESSURE 

(  .28  LB/IN) 

(  .83  LB/IN) 

(  5.1  KSI) 

.4903  N/CM 
1.4535  N/CM 
35.34  MPA 

DRAG  COEFFICIENT  OF  AIR  RESISTANCE 

1 .0000 

AXIAL  POINT  A  TO  NOSE  DISTANCE 

AXIAL  POINT  F  TO  NOSE  DISTANCE 

(  5.000  IN) 

(  2.500  IN) 

12.700  CM 
6.350  CM 

WEIGHT  OF  PROJECTILE 

(  15.598  LB) 

7.075  KG 

TRANSVERSE  MOMENT  OF  INERTIA 

(7.21  LB-INSO) 

2.110  G-MSQ 

FRICTION  COEFFICIENT  AT  BOURRELET 

FRICTION  COEFF.  AT  OBTURATOR  OR  DRIVING  BAND 

.1500 

.0400 

INITIAL  SHELL  ORIENTATION,  BETA 

THETA 

(  0.0  DEG) 

(-.1176  DEG)- 

0.000  RAD 
2.0507  MRAD 

INITIAL  C-  G.  POSITION,  X 

Y 

(  2.3105  IN) 

(  -.0147  IN) 

5  869  CM 
-.037  CM 

INITIAL  COMPUTATION  TIME  INTERVAL 

.00050  MS 

DIMENSIONS  AND  PROPERTIES  OF  GUN  TUBE 

TUBE  INCLINATION 

BORE  DIAMETER 

TRAVEL  Or  PROJECTILE 

(  0.00  DEG) 

(  4.734  IN) 

(187.000  IN) 

0.000  RAD 
12.024  CM 
4.750  M 

TABLE  3 


BASE  PRESSURE  DATA 


IN-BORE  MOTION  ANALYSIS 
EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALLOT2 


N3 

O' 


CHARGE  WEIGHT  =  .113  KGS  OR  .250  LBS 


PRESSUR 

E  ACTING  ARE<  =  1 1 3.5570 

CM. SO.  OR  17 

.6014  IN. SO. 

PRESSURE 

CENTER  TO  C.G.  DISTANCE 

=  5.8687  CM 

OR  2.3105  IN 

TIME 

FIRING  PRESSURE 

TIME 

FIRING  PRESSURE 

MILLISECONDS 

MEGAPA( PSI ) 

MILLISECONDS 

MEGAP A ( PSI ) 

0.0000 

.69 

(  100.) 

3.1000 

387.55 

(56209. ) 

.  1  coo 

1  .82 

(  264.) 

3.2000 

373.42 

(54160 . ) 

.2000 

3.05 

(  443.) 

3.3000 

358.36 

(51975. ) 

.3000 

4.66 

(  676.) 

3.4000 

342.88 

(49731 . ) 

.4000 

6.75 

(  979.) 

3.5000 

327.41 

(47487 . ) 

.5000 

9.45 

(  1371.) 

3.6000 

312.22 

(45284 . ) 

.6000 

12.94 

(  1877.) 

3.7000 

294.06 

(42650 . ) 

.7000 

17.34 

(  2515.) 

3.8000 

273.38 

(39651 . ) 

.8000 

22.99 

(  3335.) 

3.9000 

252.95 

(36687  .  ) 

.9000 

30.13 

(  4370.) 

4.0000 

233.52 

(33869. ) 

1  .0000 

39.09 

(  5670.) 

4. 1000 

215.45 

(31248. ) 

1  .  1000 

50.25 

(  7288.) 

4.2000 

198.84 

(28839. ) 

1  .2000 

64.00 

(  9233.) 

4.30C0 

183.79 

(26656 . ) 

1 . 3000 

80.45 

( 1 1669. ) 

4.4000 

170.22 

(24688. ) 

1.4160 

101 .67 

(  14746.  ) 

4.5000 

157.98 

(22913. ) 

1.4161 

101.33 

( 1 4697 . ) 

4.6000 

146.93 

(21310. ) 

1 .5000 

119.05 

(17267. ) 

4.7000 

136.94 

(19861 . j 

1 . 6000 

143.09 

(20754.  ) 

4.6000 

127.89 

(  18549.  ) 

1 .7000 

170.78 

(24770.  ) 

4.9000 

119.67 

( 17357  .  ) 

1 .8000 

202.56 

(29379. ) 

5.0000 

112.21 

( 16274  .  ) 

1 .9000 

236.44 

(34293.  ) 

5. 1000 

105.40 

(15287.) 

2. 0000 

270.36 

(39285.  ) 

5.1050 

105.64 

(  15322.  ) 

2 . 1 000 

305.45 

(44302.  ) 

5.2000 

99.30 

(14402.) 

2.2000 

337.66 

(40974. ) 

5.3000 

93.53 

(13565.) 

2.3000 

365.79 

(53053.  ) 

5.4000 

88.27 

(12802. ) 

2.4000 

388.43 

(56337. ) 

5.5000 

83.45 

(12104.) 

2.5000 

404.75 

(58704. ) 

5.6000 

79.03 

(11463.) 

2.6000 

414.50 

(60118.) 

5.7000 

74.97 

( 10874. ) 

2.7000 

418.03 

(60630. ) 

5.8000 

71.23 

(10331 . ) 

2.8000 

416.08 

(60348.  ) 

5.9000 

67.76 

(  9828.) 

2.9000 

3.0000 

409.66 

399.71 

(59416. ) 
(57973.  ) 

5.9450 

66.29 

(  9614.) 

TABLE  4 


OBTURATOR-TUBE  PRESSURE  DATA 


IN-BORE  MOTION  ANALYSIS. 

EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALLOT2 

FRICTION  COEFFICIENT  USED  IN  FRICTION  FORCE  COMPUTATION, 
AT  OBTURATOR  =  .04000 

AT  BOURRELET  =  .15000 

OBTURATOR-TUBE  PRESSURE  VS.  MUZZLE  DISTANCE 


MUZZLE  DISTANCE  PRESSURE  MUZZLE  DISTANCE  PRESSURE 

METERS  (INCHES)  MEGAPA(KSI)  METERS  (INCHES)  MEGAPA(KSI) 


0.000  (  0.00)  0.00  (  0.00)  4.750  (187.00)  70.67  (  10.25) 

AVERAGE  OBTURATOR-TUBE  PRESSURE  =  35.34  MEGAPA  OR  5.13  KSI 


TABLE  5 


*****SOM£  PRELIMINARILY  COMPUTED  VALUES  OR  GIVEN  DATA***** 


GIVEN  AND  COMPUTED  C.G.  POSITION  WITH  RESPECT  TO  ROTATING  BAND  OR  OBTURATOR  CENTER 
X  =  5.&6361  CM,  OR  2.31047  IN. 

Y  *  -.03743  CM,  OR  -.01474  IN. 

GIVEN  INITIAL  THETA  =  -.0020507  RAD.  OR  -.1174961  DEG. 

COMPUTED  MIN.  CONTACT  ANGLE  =  -.  0020507  RAD ,  OR  -.1174961  DEG 


C.G.  ECCENTRICITY  (DISTANCE  TO  PROJ.  AXIS)  «  -.02540  CM,  OR  -.01000  IN. 

C.G.  ECCENTRICITY  (UNBALANCED  MOMENT)  =  -1.76222  CM-KG,  OR  -2.49563  IN-OZ 


INITIAL  OR  CONSTANT  TIME  STEP  SIZE  =  .000500  MS 

TIME  STEP  SIZE  INCREMENT  FOR  EACH  COMPUTATION  =0.000000  MS 

EVERY  20TH  COMPUTED  DATA  ARE  RECORDED  IF  WANTED 


N> 

00 


NC 

1 

301 

601 

901 

1201 

1501 

1801 

2101 

2401 

2701 

3001 

3301 

3601 

3901 

4201 


TABLE  6 


*****  INTERMEDIATE  PRINTOUT  of  computed  values  of  main  quantities  vs. 


TIME***** 


(DATA  IN  ENGLISH  UNITS  ARE  PRINTED  OUT  ON  NEXT  LINE  WITH  E12.5  FORMAT) 


' IME , MS  X.CM(IN) 

THETA.MRD  Y,CM(IN) 

YAW  V ,  RPS 

XDT,MPS(FPS> 

XDDT , GEE 

YDT.MPS(FPS) 

YOOT.GEE  ON  ,  N  ( l  B ) 

YAW  A.RPSS  BN.N(LB) 

0.0000  0.000 
0. 

-2.051  -.097 

14738E-01 

0.000 

0,00 

0. 

0.00 

0.00 

0. 

0.00  0. 

0. 

0.  1.1B-00 

.246696-08 

.1500  .000 
•89383E-05, 

-2.050  -.037 

- .  1 4739E-01 

.037 

.02 

. 54536E-0 1 

82.26 

-.00 
- . 22869E-02 

-1.82  0. 

0. 

5.65E+02  0. 

0. 

.3000  .002 
.94453E-03 

-2.053  -.037 

-.  14751E-01 

-.168 

.39 

. 12632E+01 

446.13 

-.00 
-.1 1 206E-01 

-2.12  9.0E+00 

.2031 1E+01 

-2.22E+03  2.2E+01 

.493196+01 

.4500  .015 

.59455E-02 

-2.093  -.038 

- . 1 4782E-01 

-.225 

1 .44 
. 471 72E+0 1 

1008.81 

-.01 
- .  247S5E-01 

-3.81  -1.8E+02 
- . 39855E+02 

2.16E+03  1.7E+02 

.393056+02 

.6000  .050 

. 19865E-01 

-2.112  -.038 

-. 14842E-01 

-.187 

3.47 
.  1  1398E+02 

1800.76 

-.01 
-.41 698E-01 

-2.80  -1 . 2E+02 
- . 2  7496E+02 

-2.88E+03  3.1E+02 

.69151 E+02 

.7500  .126 

.49789E-01 

-2.175  -.038 

-. 14930E-01 

-.550 

6.94 
. 22783E+02 

2982.68 

-.02 
- . 57253E-01 

-4.32  -3.3E+02 
-.741556+02 

3.16E+02  6.0  E+.02 

.  13430E-HJ3 

.9000  .269 

. 10593E+00 

-2.243  -.038 

-. 15054E-01 

-.407 

12.47 
.  40924E  +  02 

4612.56 

02 

-.81 200E-01 

-4.90  -4.5E+02 
-.1015*6+03 

-1.17E+Q3  9.5E+02 

.21 329E+03 

1 .0500  .515 

.  20235E+00 

-2.341  -.039 

-. 15223E-01 

-.930 

20.92 
.  68650E+02 

6991  44 

-.03 

-.107136+00 

-6.89  -6.6E+02 
-.  14864E+03 

-2.73E+03  1  . 4  E  +  03 

.32362E+03 

1.2000  .917 

. 361 04E+00 

-2.486  -.039 

-. 15453E-01 

-.933 

33.44 
.  10970E+03 

10154.09 

-.05 
-.1 52226+00 

-11.22  -1 . 0E+03 
- .  23551 E+03 

1.59E+02  2.1 E+03 

. 48220  E+03 

1.3500  1.546 

.60848E-1-00 

-2.656  -.040 

-. 15775E-01 

-1 .468 

51.38 
. 1 685QE+03 

1 *34 1 .44 

-.07 
-.21 3386+00 

-14.65  -1 . 3E+03 
- . 2  8944  E+03 

-4.78E+03  3. 0E+03 

.  68384  E+03 

1.5000  2.491 

. 98084E+00 

-2.901  -.041 

1 6235E-01 

- 1 . 655 

75.90 
.  24902 E+03 

19159.22 

-.  09 
- . 297486+00 

-19.63  -1 . 7E+  03 
-.379406+03 

-4.18E+01  4 . 3  E+03 

.  97023E  +  03 

1.6500  3.863 

. 1 5209E+C 1 

-3.187  -.043 

16860E-01 

-2.3^4 

108.50 
.  35598E  +  03 

25356.88 

12 

- . 399726+00 

-23.61  - 1 .86+03 
- .  40658E+03 

-7.20E+03  5 . 9  E  +  03 

.  13302E+04 

1 .6000  5 . 797 

. 22821 £+01 

-3.592  -.045 

-. 17697E-01 

-2.391 

151.14 
. 49585E+03 

32B19.09 

16 

-.538506+00 

-34.15  -2 . 26+03 
-  . 4872 1 E+03 

-9 . 62E+02  8 . 1  E  +  03 

.  18268E  +  04 

1 . 9500  8 .456 

. 33292E+01 

-4.062  -.048 

-. 18826E-01 

-3.571 

205 . 55 
.  67436E+03 

41174.94 

-.22 

-.721896+00 

-41.36  -2  .  OE+03 
- .  44763E+03 

-7.83E+03  1.1E+04 

•24472E+04 

2.1000  12.025 

. 47342E+0 1 

-4.670  -.052 

- . 2031 7E-01 

-4. 380 

272.34 
.89351 E+03 

49644.62 

-.29 

-.941496+00 

-49 .29-1 . 7E+03 
-.371506+03 

-2.69E+03  1.4E+04 

.325826+04 

45C 1 

2.2500  16.687 

.656956+01 

-5.365  -.056 

- .  22225C-0 1 

-5.036 

4801 

2.4000  22.r,08 

. 8S00UE  +  0 1 

-6.156  -.062 

-.24525E-01 

-5.204 

5101 

2.550C  29.918 

.  1  1 7796+02 

-6.993  -.039 

- .  27087E-01 

-5.995 

5401 

2.7000  38.694 

. 15234E+02 

-7.803  -.075 

-  .  2961 8E-0 1 

-4.524 

5701 

2.8500  48  966 

. 19278E+02 

-8.424  -.080 

-.31 68eE-01 

-3.598 

6001 

3.0000  50.717 

.239046+02 

-8.714  -.083 

-.32787E-01 

146 

6301 

3.1500  73.Q99 

. 29094E  *02 

-8.555  -.082 

-.32464E-01 

2.402 

6601 

3.3000  88.443 

. 34820E+02 

-7.869  -.077 

- . 30488E-0 1 

6.579 

6901 

3.4500  104.272 

. 4 1 052E+02 

-6.741  -.069 

- . 26982E-0 1 

8-379 

7201 

3.6000  121.301 

.47756E+02 

-5.314  -.057 

- . 22495E-01 

10.219 

7501 

3.7500  139.445 

. 54900E+02 

-3.904  -.045 

-. 17904E-01 

8.578 

7801 

3 .9000  158.604 

. 624426+02 

-2.693  -.036 

-. 14213E-01 

6.809 

8101 

4.0500  178.668 

.  703426+02 

-1.908  -.031 

-. 12190E-01 

4.996 

8401 

4.2000  199.535 

. 78557E+02 

-1.442  -.030 

-. 1 1765E-01 

-.527 

8701 

4.3500  221.112 

. 87052E+02 

-1.938  -.032 

-. 12565E-01 

-4.073 

9001 

4.5000  243.321 

. 95796E^02 

-2.562  -.037 

-. 14679E-01 

-5.576 

9301 

4.6500  266.092 

.  1 0476E  +  03 

-3.430  -.045 

-. 17713E-01 

-4.960 

9601 

4.8000  289.368 

.  1  1392E  +  03 

-4.087  -.052 

- . 20569E-01 

-4.358 

9901 

4,9500  313.097 

,123276+03 

-4.564  -.056 

-.22086E-01 

-.739 

10201 

5.1000  337.237 

. 132776+03 

-4.294  -.055 

- . 21 648E-01 

3.242 

10501 

5.2500  361.749 

-3.700  -.049 

5.447 

351.10 
1 15196+04 

439.87 
14431 6+04 

535.63 
1 75736+04 

634.90 

20830E+04 

734.46 

24096E+04 

831.75 
. 272886+04 

925.07 
.  303506+04 

1013.36 
. 332476+04 

1096.17 
.35963E« 04 

1173.43 

.384986+04 

1244.64 

.408356+04 

1308.63 
.  42934E+04 

1365.45 

.447986+04 

1415.76 

.464496+04 

1460.35 
.479126+04 

1500.04 

.492146+04 

1535.51 

.503786+04 

1567.36 
.514236+04 

1596.06 

.523656+04 

1622.10 
*  532166+04 

1645.85 


57207.Se 

-.36 
-.1 17476+01 

-45.17  -5.16+02 
-.115526+03 

-6.77E+OJ  1.9E+04 

.42420E+04 

63200.56 

-.42 

-.137336+01 

-32.98  9.26+02 

.206706+03 

-9.96E+02  2.4E.04 

•53989E.04 

G6651 .54 

-.44 

-.144486+01 

2.57  3.26+03 

,712186+03 

8.39E+02  3.0E+04 

. 66586E+04 

68008.11 

-.40 
- . 1 3279E+01 

.  51.57  5.66+03 

.125706+04 

7.83E+03  3.SE+04 

.70881 E*04 

67146.43 

-.28 

-.923126+00 

112.30  7.96f03 

.177536+04 

1.51E+04  3.9E+04 

. 88694E+04 

64979.08 

-.08 
- . 25451 E+00 

165.83  9.5E+03 

.2 1 2766+04 

1.94E+04  4.2E+04 

. 93662 E* 04 

61826.45 

.19 

.635726+00 

193.44  9.66+03 

.216836+04 

2 . 48E+04  4.1 £+04 

. 91 756E+04 

58207.54 

.  47 

.154916+01 

180.89  8. 36+03 

.105636+04 

1.89E.04  3.6E+04 

.  81 826E+04 

54420.47 

.  70 

.229416+01 

116.31  5.06+03 

.1 1 268E+04 

1 . 49E+04  2.9E+04 

. 84742 E+04 

50685.55 

.80 

.260866+01 

12.58  6.96+02 

.155896+03 

-1.36E+03  1.9E+04 

•42955E+04 

46035.12 

.73 

.239526+01 

-106.29  -4.06+03 
-.907376+03 

-4.56Ef03  9.3E+03 

. 20957 E+04 

41010.20 

.50 

.162926+01 

-201 .39  -7.46+03 
-.165956+04 

-3.11E+04  1.2E+03 

. 27537 E+03 

36351 .68 

.  19 

.633526+00 

-189.14  -8.26+03 
- . 1 8398E+04 

7.09E+02  0. 

0. 

32155.63 

-.03 
-.1 10986+00 

-131 .14  -5.86+03 
-.129826+04 

-5.93E+04  0. 

0. 

28562.59 

-.25 

-.804876+00 

-150.58  -7.06+03 
- . 1 5Q32E+04 

1.06E.04  0. 

0. 

25470.63 

-.46 

-.151026+01 

-1 12.97  -4.6E+03 
-.102566+04 

-2.19E+04  1.4E+03 

.  31128E+03 

22837.14 

-.53 

-.175516+01 

19.15  -7.46+02 
-.167476+03 

1.54E.04  7.7E+03 

. 17236E.04 

20530.91 

-.40 
-. 13014E+01 

161.42  4.1E+03 

.916786+03 

"8.71 E+02  1.3E+04 

2850. £.04 

18573.40 

10 

-.330566+00 

229.43  5.66+03 

.125806+04 

4.36E+04  1.6E.04 

. 36069E.04 

1 6858.74 

.25 

.821626+00 

233.23  6.36+03 

•  1 41 686+04 

8.02I+0*  1.SC-.04 

.33024C.04 

15408.53 

.53 

128.90  2.66+03 

3.03E+04  1.0E+04 

.  1 4242E+03 

- . 1 9275E-01 

. 53998  E+04 

.1 7295E+01 

. 591 90E+03 

. 23456E+04 

10801 

5.4000  386.602 

.  15221  E+03 

-2.630  -.040 

-.  15802E-0I 

7.512 

1667.53 
. 547096+04 

14101 .60 

.62 

.201 87E+0 1 

-9.45  -9.7E+02 
-.21 829E+03 

-1.61 E+04  3 . 0  E+03 

.67997 E+03 

11101 

5.5500  411.766 

.16211 E+03 

-1.639  -.032 

-.  1  2402E-0 1 

7.969 

1637.42 
.55361 E+04 

1 2965.50 

.  51 

. 1 689 1 E+0 1 

-96.32  -5 . 1 E+03 
-.1 1 460E+04 

5.36E+04  0. 

0. 

11401 

5.7000  437.217 

. 17213E+03 

-.044  -.025 

- . 97092E-02 

10.254 

1705.72 
. 55962  E+04 

11953.00 

.  42 

. 1 3681 E+0 1 

-22.64  -1 . 5E+03 
- .  32825E+03 

-5.35E+04  0. 

0. 

11701 

5.8500  462.931 

. 13226Ef03 

.798  -.019 

- . 73029E-02 

2.678 

1722.64 
. 5651 7E+04 

1 1066.80 

.  39 

. 1 2908E+0 1 

-37.13  -3. 1E+03 
- .  7 0472E+03 

6.91 E+03  0. 

0. 

11841 

5.9200  475.016 

. 18701 E+03 

1.05S  -.016 

- . 62627E-02 

5.445 

1730. 10 
.  56762E+04 

1 0680.20 

.36 

.1 1 777E+01 

-57.92  -4.8E+03 

1 0688E+04 

6 . 1 0E+0 0. 

TOTAL  NO  OF  COMPUTATION  =•  11841 


TOTAL  NO  OF  DATA  RECORDED  ON  TAPE  -  593 


TABLE  9 


»****D£T AI LED  OUTPUT  OF  COMPUTATION***** 

BASE  PRESSURE,  NORMAL  AND  FRICTIONAL  FORCES  AND  AIR  RESISTANCE 


BP  -  BASE  PRESSURE,  KPA  (.145  PSI) 

BN  =  BCURRELET  NORMAL  FORCE,  N  (LB) 

FD  =  AERODYNAMIC  DRAG  FORCE,  N  (LB) 

FL  -  AERODYNAMIC  LIFT  FORCE,  N  (LB) 

ON  =  BAND  OR  OBTURATOR  NORMAL  FORCE,  N  (LB) 

FOI  =  BAND  OR  OBTURATOR  FRICTIONAL  FORCE  FROM  BAND-TUBE  PR.,  N  (LB) 
FQ2  =  BAND  OR  OBTURATOR  FRICTIONAL  FORCE  FROM  ON,  N  (LB) 

FB  *  BCURRELET  FRICTIONAL  FORCE,  N  (LB) 

FX,  FY  =  X,  Y  TOTAL  FORCE  AT  C.G.,  N  (LB) 

1  N  =  .22482  LB 


uj 


NO 

TIME, MS 

BP, KPA 

ON  ,  N 

8N  ,  N 

1 

0 . 0000 

6 . 395E+02 

0  . 

1  .  106E-08 

0. 

21 

.  2000 

3 . 054E+03 

1 . 731CE+01 

0. 

0. 

41 

.4000 

6 . 750E+03 

-1  .1011 E  +  02 

1 . 193E+02 

o. 

61 

.6000 

1 .2946+04 

-1  . 2230  E  +  02 

3 . 076E+02 

0. 

81 

.  6000 

2.2S9E+04 

-4 . 1 860E+02 

7. 187E+02 

o. 

101 

1  . 0000 

3. 909E+04 

-5.2775E+02 

1 . 24 1 E+03 

o. 

121 

1  .2000 

6 . 400E  +  0  4 

-1  . 0475  E  +  03 

2.  145E  +  03 

o. 

141 

1 .4000 

9 . 874E+04 

~ 1  . 4546  E  +  03 

3 . 436E+03 

o. 

161 

1 .6000 

1 .431 E+05 

-1  . 7334  E  +  03 

5.320E+03 

o. 

181 

1 .8000 

2.026E+0S 

-2. 1671 E+03 

8. 126E+03 

o. 

201 

2 . 0000 

2.709E+05 

-1  . 9 1 48  E  +  03 

1  .  199E+04 

0. 

221 

2.2000 

3.377E+05 

“9 . 6670  E  +  02 

1  . 732E+04 

o. 

241 

2.4000 

3.S84E+05 

9 . 1 942  E  +  02 

2.401 E+04 

0. 

261 

2.6000 

4. 145E+05 

3 . 8245  E  +  03 

3. 154E+04 

0. 

281 

2.8000 

4 . 1 6 1 E  +  05 

7 . 3  305  E  +  03 

3.814E+04 

0. 

301 

3.0000 

3 . 997E+05 

9 . 4  635  E  +  03 

4. 166E+04 

o. 

321 

3.2000 

3.734E+05 

9 . 2228  E  +  03 

3.97SE+04 

0. 

341 

3.4000 

3.429E+05 

6 . 4082  E  +  03 

3.15SE+04 

0- 

361 

3.6000 

3 . 1 22E+C5 

6 . 9342E+02 

1.911 E+04 

0. 

381 

3-8000 

2 . 734E+Q5 

-5.5226E+03 

6 . 377E+03 

0. 

401 

4.0000 

2.335E+05 

-8.5068E+03 

0. 

0. 

421 

4.2000 

1  . 988E  +  0  5 

"5 . 7  74  2  E  +  03 

0. 

0. 

441 

4 . 4000 

1 . 702E+05 

~  6 . 6089 E+03 

0. 

o. 

461 

4.6000 

1 . 469E+05 

“2 . 351 7E+03 

5 . 578E  +  03 

0. 

481 

4.8000 

1 .279E+05 

4 .077 9 E+03 

1  .312E+04 

0. 

501 

5.0000 

1 . 122E+G5 

6 . 0484  E  +  03 

1  . 647  E+04 

0. 

521 

5*2000 

9.930E+04 

4 . 1 734  E  +  03 

1  . 232E+04 

o. 

541 

5.4000 

8 . 827E+04 

-9. 7095E+02 

3 . 024E+03 

0. 

561 

5.6000 

7 . 903E+04 

-4.71 30  E  +  03 

0. 

0. 

581 

5.8000 

7. 123E+04 

-1 . 4943 E+03 

0. 

0. 

593 

5.9200 

6.71 1E+04 

-4.7539E+03 

0. 

0. 

F  L , N  FD , N  FOI , N  F02 , N  F8,N  FX,N  '  FY  ,N 


0. 

-4.571E-05 
-6.861E-03 
-8.356E-02 
-5. 032E-01 
-2. 1 77E+00 
-7. 756E+00 
-2.398E+01 
-6.481E+01 
-1 . 586E+02 
-3.561 £+02 
-7. 275E+02 
-1 . 344E+03 
-2. 246E+03 
-3.420E+03 
-4.81 0E+03 
-6. 343E+03 
-7. 950E+03 
-9. 577E+03 
-1.11 6E+04 
-1 . 263E+04 
-1 . 394E+04 
-1 . 51 2E+04 
-1 . 616E+04 
-1 . 709E+04 
-1 . 792E+04 
-1 . 867E+04 
-1 . 935E+04 
-1 . 996E+04 
-2. 052E+04 
-2. 083E+04 


-2.1 969E+04 
-2. 1 969E+04 
-2 . 1 969E+04 
-2. 1967E+04 
-2.1 961 E+04 
-2. 1950E+04 
-2. 1 927E+04 
-2. 1 8B5E+04 
-2 . 1 814E+04 
-2.1701 £+04 
-2 . 1 529E+04 
-2.1276E+04 
-2 . 0924E+04 
-2 . 0459E+04 
-1  . 9872E+04 
-1  .  9 1 63E  +  04 
-1  . 3336E+04 
-1  . 7 399E+04 
-I  .6361 E+04 
-1  . 5232E  +  04 
-1  . 4022E+04 
-1  . 2743E+04 
-1 . 1 406E+04 
-1  .001 9E+04 
-8.5887E+03 
-7.121 1 E+03 
-5 . 6208E+03 
-4.091 6E+03 
-2 . 5368E+03 
-9 . 5903E+02 
-2 .3551 E+00 


0. 

-6 . 924E-0 1 
-4.405E+00 
-4 . 892E+00 
-1  .674E+01 
-2.111 E+01 
-4.1 90E+01 
“5.81 8E+01 
-6.934E+01 
“8 . 668E+0 1 
-7 .659E+01 
-3.867E+01 
-3 . 678E  +  0 1 
-1  . 530E  +  02 
-2 . 932E+02 
-3.785E+02 
-3 . 689E+02 
-2 .563E+02 
-2 . 774E+01 
-2 . 209E+02 
-3 . 403E+02 
-2.31 0E+02 
-2 . 644E+02 
-9 . 407  E  +  01 
-1 .631 E+02 
-2.41 9E+02 
-1  . 669E+02 
-3.884E+01 
-1  .885E+02 
-5 . 977E+0 1 
-1 . 902E+02 


-1  .6593E-09 

0. 

-1 .7901 E+01 
-4 . 61 38E+01 
-1 .0781 E+02 
-1 .8621 E+02 
-3.21 73E+02 
-5 . 1 546E+02 
-7 . 97 99E+02 
-1 . 21 89E+03 
-1 . 7987E+03 
-2 . 5977E+03 
-3.6021 E+03 
-4.731 1 E+03 
-5.7207E+03 
-6.2491 E+03 
-5 . 9674E+03 
-4 . 7386E+03 
-2 . 8659E  +  03 
-9.5650E+02 
0. 

0. 

0. 

-8.3668E+02 
-1 . 9683E+03 
-2.4698E+03 
-1 . 8482E+C3 
-4.5367E+02 
0. 

0. 

0. 


0. 

1 .271 E+04 
5 . 466E+04 
1 . 249E+05 
2 . 390E+05 
4.21 8E+05 
7.045E+05 
1 . 099E+06 
1  . 602E+06 
2.277E+06 
3 . 052E+06 
3.81 0E+06 
4 . 385E+06 
4 . 679E+06 
4 . 695E+06 
4.508E+06 
4 . 209E+06 
3 . 863E+06 
3.51 6E+06 
3 . 077E+06 
2 . 625E+06 
2 . 231 E+06 
1 . 906E+06 
1 . 64 1 E+06 
1 . 424E+06 
1 . 246E+06 
1 . 101 E+06 
9 . 783E+05 
8 . 748E+05 
7 . 873E+05 
7.41 0E+05 


0. 

-1 . 2309E+02 
-2. 1949E+02 
-1 .9451 E+02 
-3.4395E+02 
-3. 7591 E+02 
-7.7851 E+02 
-1  .  1527E+03 
-1 .4831 E+03 
-2.3692E+03 
-3.0584E+03 
-3.3476E+03 
-2.2878E+03 
9 . 8349E+02 
6 . 5236E+03 
1 . 1 505E+04 
1 . 3308E+04 
1 . 0100E+04 
8.7289E-+02 
-1 . 0028E+04 
“1 .4291 E+04 
“9 . 0980E+03 
-1 . 0789E+04 
-2.1 127E+03 
1  .  1 199E+04 
1 . 6649E+04 
1  . 1982E+04 
“6 . 5568E+02 
“5 . 8395E+03 
-1 .0402E+03 
-4.0181E+03 


TABLE  10 


**  *  *  *  DET  A  I  LED  OUTPUT  OF  COMPUTATION***** 

ANGULAR  DISPLACEMENTS,  VELOCITIES  AND  ACCELERATIONS  AND  MOMENTS 


THETA  =  YAW  ANGLE,  trt  I LL I R AD  I ANS 
THETAD  =  YAW  VELOCITY,  RAD/SEC 
THT ADD  =  YAW  ACCELERATION,  RAD/SEC/SEC 
ZM  =  TOTAL  MOMENT  ABOUT  fc.G.,  N-M  (IN-LB) 

BM  =  MOMENT  ABOUT  C.G.  FROM  BOURRELET  FORCES,  N-M  (IN-LB) 

OM1  =  MOMENT  ABOUT  C.G.  FROM  FORCE  F01,  N~M  (IN-LB) 

0M2  *  MOMENT  ABOUT  C.G.  FROM  OBTURATOR  (BAND)  FORCES,  N-M  (IN-LB) 
OM3  ■  LATERAL  ROTATION  RESISTING  MOMENT,  N-M  (IN-LB) 

AM  =  MOMENT  ABOUT  C.G.  FROM  AIR  RESISTANCE  FORCE,  N-M  (IN-LB) 

1  N-M  =  9.650748  IN-LB 


NG 

TIME, MS  THETA, MRAD 

THETAD, RPS 

THTAQD , RPSS 

OM1 , N-M 

0M2,N-M 

OM3,N“M  BM , N-M 

AM, N-M 

ZMt  N-M 

1 

0. 0000  -2.051E+00 

0. 

0. 

8 . 2245E  +  00 

0. 

0. 

5.526E-10 

0. 

0. 

21 

. 2000  -2.0485+00 

3 . 62835-02 

-7.8276+02 

8.2260E+00 

-1 . 059E+00 

0. 

0. 

6.0839E-09 

-1 . 652E+00 

41 

,4000  -2.0795+00 

-3  *  0376E-01 

7. 1906+02 

8.241 2E+00 

6.71 5E+00 

0. 

5 . 962E+00 

9.00G5E-07 

1 .517  E+00 

61 

.6000  -2. 1 12E+0Q 

-1 .8724E-01 

-2.880E+03 

8.281 6E  +  0  0 

7 . 458E+00 

0. 

1 . 537E+0 1 

1  .0803E-05 

-6. 078E+00 

81 

.8000  -2.201E+00 

-4.6546E-01 

2.01 1E+03 

8 . 3490E+0  0 

2. 552E+01 

0. 

3.590E+01 

6. 2422E-05 

4.243E+00 

101 

1 . 0000  -2 . 2995+00 

-7.51 19E-01 

-4.0005+03 

8 . 4529E+00 

3.2185+01 

0. 

6.201 E+01 

2.5759E-04 

-8 . 440  E+00 

121 

1 . 2000  -2 . 4865+00 

-9.3274E-01 

1 . 588E+02 

8.6061 E+0  0 

6.3065+01 

0. 

1 . 07 1 E+02 

8 . 3 182E-04 

3.351E-01 

141 

1 .4000  -2 . 735E+00 

—  1  . 64576+00 

-2. 033E+03 

8 . 8464E+00 

8 . 064  E  +  0 1 

0. 

1 . 71 6E+02 

2 . 2207E-03 

-4.291 E+00 

161 

1 . 6000  -3 . 078E+00 

-1 . 9970E+00 

-5. 956E+03 

9.21 45E+00 

1 .0565+02 

0. 

2 . 655E+02 

4.6864E-03 

-1 . 257E+01 

181 

1 . 8000  -3 . S92E+00 

-2*89076+00 

-9.6215+02 

9. 7534E+00 

1 .3205+02 

0. 

4 . 053E  +  02 

6.6629E-03 

-2 . 030  E+00 

201 

2.0000  -4.25CE+00 

-3 . 9466E+00 

-6. 759E+03 

1 .05415+01 

1 . 1 65E+02 

0. 

5.977E+02 

1 . 1343E-03 

-1  *  426E+0 1 

221 

2. 2000  -5. 122£t00 

-4.71 50  E  +  00 

-5. 939E+03 

1 . 1 639E+0 1 

5.8785+01 

0. 

8.6255+02 

-3.511 7E-02 

-1 . 253  E+01 

24  1 

2. 4000  -6. 1586+00 

-5 . 2040  E  +  00 

-9. 9626+02 

1 .3031 E+01 

-5.6495+01 

0. 

1 . 195E+03 

-1 .471 7E-01 

-2. 1 02  E+00 

261 

2.6000  -7 . 293E+0G 

-5.71 80  E+00 

9.2706+03 

1 .4521 E+01 

-2. 35 2 E+0 2 

0. 

1 . 568E+03 

-3 . 9630E-0 1 

1 . 956E+0 1 

281 

2.8000  -8.230E+00 

-4 . 0782E+00 

4 . 478E+03 

1 . 5686E+0 1 

-4.512-E  +  02 

0. 

1  . 894E  +  03 

-7 . 9282E-0 1 

9.451 5+00 

301 

3.0000  -8.7145+00 

-1  .  4569F.-01 

1 . 943E+04 

1 . 5958E+0 1 

-5.8286+02 

0. 

2 . 068E  +  03 

-1  . 2530E  +  00 

4. 1 00  E+0 1 

321 

3.2000  -8.4005+00 

3  >  8336E+00 

3.  1 12E  +  04 

1 . 4SC2E+0 1 

-5.678E+02 

0. 

1 . 975E+03 

-1 . 5357E+00 

6. 568E+01 

341 

3.4000  -7. 144E+00 

7 . 7672E+00 

9. 080E+03 

1 .2512E+01 

-3.941 E+02 

0. 

1 . 57  0E  +  03 

-1  . 3358E+00 

2.085E+01 

361 

3.6000  -5.314E+00 

1 .0219E+01 

-1 . 357E+03 

9 . 3553E+00 

-4.2576+01 

0. 

9.51 4E+02 

-5 . 7507E-01 

-2.8645+00 

301 

3.8000  -3 . 479E+00 

8.41 92E+00 

-3. 839E+03 

6 . 3986E+00 

3 . 363E+02 

0. 

3. 1815+02 

5 . 39G6E-0 1 

-8.101 E+00 

401 

4 . 0000  -2 . 1 535+00 

4 . 7714E+00 

3. 248E+03 

4.51 27  E+00 

5.  187C-502 

0. 

0. 

1 . 600GE  +  00 

6 . 854  E+00 

421 

4.2000-  -1  .  442E  +  00 

-5. 2699E-01 

-5. 9326+04 

3 . 8079E+00 

3.5235+02 

0. 

0. 

2.3552E+00 

-1  . 252  E+02 

441 

4. 4000  -2. 1285+00 

-3.6282E+00 

2.8106+03 

3 . 8007E+00 

4 . 0  30  E+02 

0. 

0. 

1 .941 7E+00 

5.931 5+00 

461 

4. 6000  -3. 1 61 E+00 

-5 . 7695E+C0 

1 . 452  6  r 04 

4 . 2370E+00 

1  .4336  +  02 

0. 

2 . 783E+02 

1  .091 2E  +  00 

3.065E+01 

431 

4.8000  -4.C87E+00 

-4 . 3584E+00 

-6.710E+02 

4 . 4  857  E+0  0 

-2 . S00E+C2 

0. 

6 . 543E+02 

2. 1954C-01 

-1 . 4 1 6E+00 

501 

5. 0000  -4.5485+00 

1 .3123Ef00 

3 . 598E+04 

4.01 06  E+00 

-3.710E+C2 

0. 

8.2075+02 

-2.6092E-01 

7 . 592  E+0 1 

521 

5.2000  -3.938E+00 

4 . 1 433E+00 

2. 01 6E+04 

2.891 5E+00 

-2.5596+02 

0. 

6. 144E+02 

3 . 9893E-0 1 

4. 2555+01 

541 

5. 4000  -2 . 630E+00 

7.51 1 7  E+00 

-1 .6096+04 

1  . 6436E+0  0 

5 . 9 1 8  E  +  0 1 

0. 

1 .51 0E+02 

1 . 905GE+00 

-3 . 395  E  +  0 1 

561 

5. 6000  -1.1 735+00 

1 . 0607E+01 

4.41 9E+04 

7 . 3737E-0 1 

2.8776+02 

0. 

0. 

3.681 6E+00 

9. 326E+01 

581 

5,8000  6.490E-01 

3.71 89E+00 

-4. 718E+04 

1 .9731 E-01 

9. 1 36  E  +  0 1 

0. 

0. 

5 . 9952E  +  00 

-9 . 957  E+01 

593 

5.9200  1 . 059E+00 

5 . 4448E+00 

6. 101 E+04 

3.7507E-04 

2.908E+02 

0. 

0. 

6 . 5894E+00 

1 . 288  E+02 

9 1 1 


FIGURE  1 


C.G.  MOTION,  AXIAL 


example  of  computation  of  program  BALLOT 2 

CG  ECC=-.025MM,CG  DI ST-23 . 65CM ,  OBTR-BRT  L*11.77CM 
BORE  DIAMETER  =  12.024CM,  BOURRELET  0.=  11 . 976CM 
FRICTION  COEFF.=  .040(OBTURATQR) f  . 1 50 ( BOURRELET ) 
SPRING  C.  ( MN/CM ) *  1 . 45 (OBTURATOR ) ,  . 49 ( BOURRE LET ) 


6.4000  1 - 1 - 1- 

1 

1 

1 

1 

5.6000  - 
1 
1 
1 
1 

4.8000  - 
1 
1 
1 
1 

4.0000  - 
1 
1 


■1- 


-1* 


-1- 


cccccc 
cc 
cc 
cc 
cc 
cc 


cc 

cc 

cc 

cc 

cc 

cc 

c 


1 

3.2000  - 
1 
1 
1 
1 

2.4000  - 
1 
1 
1 
1 

1 .6000  - 
1 
1 
1 
1 

.8000  - 
1 
1 

1  CCC  BBBBB  AAAAA 

1  CCCCCCBBBBBBBB AAA AAAAA 

-.0000  CCCCCCCBBBEBBBAAAAAAAA - 

0.0000  2.0000 


CCC 
CCC 
CCC 
CCC 
CCC 
CCC 
CCC  bbb 
CCBB 
BCC 
BBBCC 
BBB  CCC 
BBB 


CC  BBB 

CC  BBB 

CC  BBB 

CC  BBB 

CC  B8B 

CC  BBB  AAA 

CC  BBB  AAA 

CCC  BBB  AAA 

CC  BBB  AAA 

CC  BBB  AAA 

CC  BBB  AAA 

CC  BBB  AAAA 

CCC  BBB  AAAA 


BBBBQBB 
BBBBBB  AAA 
BBBBB  AAA 
BBBBB  AAA 

BBB8  AAA 

AAA 
AAA 
AAA 
AAA 
AAA 
CC  AAA 
CC  AAA 


- 1 

1 

1 

1 

1 

1 

1 

1 

1 

A- 
AA1 
AA  1 
BBBBB 1 

1 


CCA  AA 
CC 

AACC 
AAA  CC 


CCC 

CCC 

CCC 

CCCC 

CCCCC 


4 1 0000 


1.0000  3.0000 

TIME,  MILLISECONDS 


5.0C00 


1 
1 
1 
1 

1 
1 
1 
1 

1 
1 
1 
1 

1 
1 
1 
1 

CCCCC- 
CC 1 
1 
1 
1 

- 1 

6. 


CURVE  A  =  TRAVEL,  B  =  VELOCITY,  C  =  ACCELERATION 

SCALE  OF  VERTICAL  AXIS - 

TRAVEL, M.  VELOCITY, X400  M/SEC  ACCELERAT ION , XI 5000  GEE 


0000 


FIGURE  2 


C.G.  MOTION,  AXIAL 


U> 


EXAMPLE  OF  COMPUTATION  OF  PROGRAM  CAL LOT2 

CG  ECC  =  -  .  025MM  »  CG  DI  ST  =  28 . 6.5CM ,  QBTR“EJRT  L*11.77CM 
BORE  DIAMETER  -  12.024CM,  BOURRELET  D.=  1 1 . 976CM 
FRICTION  COEFF  .  =*  .  040  ( OBTURATOR )  f  .  1  50  (  BOURRELET  ) 
SPRING  C.(MN/CM)=1 . 45 (OBTURATOR ) #  . 49 ( BOURRE LET ) 


6,4000  1- 

1 
1 
1 
1 

5.6000  “ 

1 
1 
1 
1 

4.8000  - 
1 
1 

1 

1  cc 
4.0000  -  CC 
1  c 
1  c 

ICC 

ICC 

3.2000  -CC 
ICC 
1C 
1C 


- ! - 1 

1 

1 

AAAA 

BBBBB8BB 

BDBBBBBBBBBB 


C 

CCCCC 
CC  CCC 


CCC 

CCC 


ABBB 

AB8B 

AAABB 

AAABB 


2.4000 


1 . 6000 


.8000 


1C  AAABB 
“C AAA  3B 
1 CAA  SB 
CCA  6B 
CC  SB 
CC  83 
CC  BB 
CC  BB 
CCBB 
C  6 
CB3 
CBB 
C3 
CB 
CB 
C 


SBBBBBBBBAAA 
BBBBBBBB  AAAAA 
BBBBBB  AAAAAA 
beBbbb  AAAAAA 
BBBBBB  AAAAA 
BBBB  AAAAA 

BBBB  AAAAA 
BBBB  AAAAA 
BBBB  AAAAA 
BBB  AAAA 
CCC  BBBAAAAA 
CCCAAAA 
BBCCCA 
BBBACCC 
CC 
CCC 
CCC 
CCC 
CCC 
CCCC 
CCC 
CCCC 
CCCC 
CCCC 
CCCCC 
CCCCC 

CCCCCC 

CCCCCCCC 

CLCCCCCCC 


.0000  c- 


*  1  * 


-1- 


*r 


1 
1 
1 
1 

1 
1 
1 
1 

1 
1 
1 
1 

1 
1 
1 
1 

1 
1 
1 
1 

1 
1 
1 
1 

ccccccccccc 

cccc 

1 

1 

1 

-“-1 


-1- 


0.0000 


.8000 


1.6000  3.2000 

2.4000 

TRAVEL,  METERS 


4.8000 


4.0000 


CURVE  A  =  TIME,  3  =  VELOCITY,  C  =  ACCELERATION 

SCALE  OF  VERTICAL  AXIS - 

TIME, MS.  VELOCITY, X3CC  M/SEC  ACCE LERAT ION , XI 5000  GEE 


138 


.6000 


.4000 


.2000 


0.0000 


-.2000 


-.4000 


-.6000 


-.8000 


-1 .0000 


FIGURE  3  C.G.  MOTION,  NORMAL 

EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALLO  T2 


■  #  / 


CG  ECC=- . 025MM , CG  DI 5T=28 . 65CM,  0BTR-8RT  L 
80RE  DIAMETER  =  12.024CM,  BOURRELET  D.=  1 1 . 976CM 
FRICTIGN  COEFF . =  . 040 (OBTURATOR) f  . 1 50 ( BOURRELET ) 
SPRING  C. (MN/CM) = 1 . 45 ( OBTURATOR ) t  . 49 ( BOURRELET ) 


‘1- 


■1* 


■1- 


“1* 


■1- 


BBB 
BB  BB 
BB  BB 


1  BB  BB 

1  B  B 

CCCC  BB 

1  CCCBCC  BB 

1  CCCBB  CC  BB 

1  CC  BB  CC  BB 

1  CC  B  CC  BB 

CCCCCCCCCCCCCCCCCC  CCC  BB  CC  BB 

1  BBBBCCCCCCCCCCCCC  B  CC  BB  CC 

1  BBBB  BB  CC  BB  CC 

1  BBB  BB  CC  CCCCC 

1  BBQB  BBB  CCCCCCCC 

BBBB36B  CCCC  B3 

1  BBBB  CC  BEBB 

1  AAA  BBB 

1  AAAAAA 

AA AAA AAA AAAA  AAA  AA 

AAAAAAAAAAAAAAAAA 


BBB 
BBBBB 
CCCCB  BBB 
CCCCC 
CC  BCC 
CC  BCC 
C  BB  C 
CC  BB  CC 
CC  B  C 


BB 

BB 

BB 

BB 

B 


•1 
1 
1 
1 
1 

1 
1 
1 
1 

BBBB- 
B 1 
1 
1 
1 

CC  - 
CCCCCC1 
CCCC  Cl 
CC  1 
AA  1 
AAA¬ 
AAA  1 


aaaaa 

AAAA 

AAA 

AAA 

AA 

AA 

AA 

AA 

AAA 


AA 

AA 

AA 

AA 

AA 

AA 

AA 

AA 

AA 

AA 


AA 
AA 

AAA  AA 
AAA  AAA 
AAAAA 


AA 

AA 

AA 

AAA 

AA 


AAAAAAA 

AAAA 


-1- 


■1- 


“1- 


'.0000 


1 .0000 


2.0000  4.0000 

3.0000 

TIME,  MILLISECONDS 


5.0000 


6.0000 


CURVE  A  =  DISPLACEMENT,  B  =  VELOCITY,  C  =  ACCELERATION 
SCALE  OF  VERTICAL  AXIS - 

D I S P L  f  MM •  VELOC I T Y , X  2  M/SEC  ACCE LERAT ION , X  1000  GEE 


FIGURE  4  C.G.  MOTION,  NORMAL 

EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALL0T2 


CG  ECC=- . 025MM , CG  DI ST  =  28 . 65CM ,  OBTR-SRT  L=I1.77CM 
BORE  DIAMETER  =  12.024CM,  BOURRELET  D.«  1 1 . 976CM 
FRICTION  COEFF.r-  . 040 (OBTURATOR) ,  .  1  50  (  BOURRELET) 
SPRING  C. ( MN/CM ) = 1 . 45 ( OBTURATOR) f  .49 ( BOURRELET) 


.6000 


M* 


•1* 


—  1* 


■1- 


-1- 


.4000 


» 2000  - 
1 
1 
1 
1 


BBBB 
BB  BBB 
SB 
BB 
BB 
BB 


uo 

\o 


BB 
BB 
BB 

ccccc 

CC  BCCC  BB 

CC  BB  CCC  BE 

CC  BB  CC  BB 

CC  BB  CC  BB 

0.0000  CC  CC  B  CC  BB  CC 

CCCC  BB  CC  BS  CCC 

BB  BB  CCC  BB  CCC 

IBB  BB  CC  CCCCBCCC 

IBB  BB  CC  CCCCCCCCC 

BB 


BBBBB 

BEBBBBBB 


CCCCCCBB 
CCCC  CCC 


BBBB 


CCC 

CC 

CC 

CC 


BBCC 
BB  CC 
BB 

BB 
B 


CC 

CCC 

CC 


BBBBBBBB 

BB 

1 

1 

1 


-.2000 


-  BBBBB 
1  BBB 
1 


BS 

BB 

BBB 

BBB 

BBB 


*.4000 


A 

AA 
AA 
1  AA 
1  AA 


1 

-.6000  - 
1 
1 
1 
1 

-.8000  - 
1 
1 
1 
1 

-1 .0000  1- 


AA 
AA 
AAA 
AA 

AA  AA 

AA  AA 

AA  AA 

AA  AAA 

AA  AA 

AAA  AA 
AAA  AAA 
AAAAA 


CCCCC 
CC  BBBBBBB 

AAAAA  BBBB 
A A A AAA A AAA 
AAA  AAAA 

AAA 
AAA 
AAAA 
AAAA 


CCCC  - 
CCC  CCC  CCCC 
CCCCC C  CC 

CC  1 

AAA 
AAAA- 
AAAA  1 


AAA 

AAA 

AAA 


AAA 
AAA 
AAAA 
AAAA 
AAAA AAA 


“  1 ' 


■1- 


-1- 


o.oooc 


,8000 


1 .6000  3.2000 

2.4000 

TRAVEL,  METERS 


4.8000 


4.0000 


CURVE  A  = 


DISPLACEMENT, 
SCALE  0? 
VELOCI TY,X 


a  =  VELOCITY,  C  =  ACCELERATION 

VERTICAL  AXIS - 

2  M/SEC  ACCELERATION^  1000  GEE 


d:spl,mm. 


FIGURE  5  YAW  ANGLE  AND  VELOCITY 

EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALLQT2 


O 


CG  ECC=- . 025MM , CG  DlST^28 . 65CM,  08TR-BRT  L=11.77CM 
30RE  DIAMETER  =  12.024CM,  BCURRELET  D.=  11.976CM 
FRICTION  COEFF.=  . 040 ( GB  TURATOR ) t  . 1 50 ( BOURRELET ) 
SPRING  C. (MN/CM) = 1 . 45 { 03TURAT0R ) f  . 49 ( BO JRRELET ) 


6 . 0000  1 - 1 - 1  — 

1 

1 

1 

1 

1 
1 
1 
1 

1 

1 

1 

1 

O.COOO  BBBBeBBBBBB 

1  BB  BBB9BBBB 
1  BBBB 

1  BBB 

1  BBBB 

-2.0000  AAAAAAAAAAA  BBBB 


4.0000 


2.0000 


BBB 
BBB 
BB  BB 


BB 

BB 

B 


BB 

BB 

B 

B 

BB 

B 


BB 

B 

B 

BB 

BBB 

B 


- BB¬ 
BB 
BB 
B3 
BBB 

BBBBBB 
BBBBBB 
BBB  BB 
BB  B 


BB 

BB 


1 
1 

BB  1 
BB1 
BB¬ 


BB 

BB 

BB 

BB 

B 


-4.0000 


“6.0000 


1 
1 
1 
1 

1 
1 
1 
1 

1 

1 

1 

1 

-8.0000  - 
1 
1 
1 
1 

1- 


AAAAAAAA 

AAAAA 


AAAA 

AAA 

AAA 

AA 


BB 
8BB 
BBBBB88B 
BBBBBB 
BB 


AA 

AA 

AA 

AA 

A 


BB 

BBB 

BB  1 

BB 

BB 

BB  1 

B  3 

B 

AA  1 

B 

B 

AA  1 

B 

BB 

A  - 

B 

BB 

AA  1 

B 

BB 

AA  1 

BB 

B 

A  1 

AABBB 

B  AA  1 

AA  BBB  BBB  AA 

- 

iA  BB3BB3 

AA 

1 

i  B6B 

AA 

1 

i  BBA 

AA 

1 

AAA 

AA 

1 

AA 

aaa 

AAA 

AA 

AA 

AA 

AA 

AA 

AA 


AA 

A 

AA 

AA 

AA 

AA 

AA 

AA 

AA 


AAA  AAA 
AAAAA 
AA 


AAAAAA 

AAAA 


-10.0000 


•1' 


0.0000 


1 .0000 


2.0000 


TIME, 


3.0000 
MILLISECONDS 


—  1 - 

4.0000 


5.0000 


6.0000 


CURVE  A  = 
YAW  ANGLE 


YAW  ANGLE,  B  ==  YAW  VELOCITY 

SCALE  OF  VERTICAL  AXIS - 

MILLIRAD  YAW  VELOCITY,  X  2  RAD/SEC 


FIGURE  6  YAW  ANGLE  AND  VELOCITY 

EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALLOT2 


6.0000 


4 . COCO 


2.0000 


0.0000 


-2.0000 


-4.0000 


-6. 0000 


-8.0000 


-10.0000 

0 


CG  ECCs-.02SMM.CG  DI ST  =  28 . 65CM ,  OBTR-BRT  1-11. 77CM 
BORE  DIAMETER  =  12.024CM,  BOURRELET  D.=  i 1 . 976CM 
FRICTION  COEFF . =  • 040 (OBTURATOR ) ,  . 1 50 ( BOURRE LET ) 
SPRING  C. (MN/CM)*1 . 45 ( OBTURATOR ) f  . 49 ( 30URRE LET ) 


1 — 

1 
1 
1 

1 

1 
1 
1 
1 

1 
1 
1 
1 

B 
B 

BB 
BB 
BB 
ABB 

abbbbb 
aabbbb 
AA  BB 
AA 
-A 
1  AA 
1  AA 
1  A 
1  AA 
-  AA 
1  AA 
1  AA 
1  A 


-1- 


BBB 
B68B 
B0  6808 
BB  BB 
BB  B8 

B  B 

BB  BB 

B  BBBB 

B  BB 

BB  B 

bb  b 

BB  BB 

BB  BB 

BB  B 

BB  B 

9  B 

ABB 

AAAABABB 
AAA  BBSBB 
AAA 
AA 


—  BB— 
BBB 
BBBB 
B  B 


B 

BB  BB 
BBBB  B 
BB  BBB 


B 

BB 

BB 

B 


BB 

BB 

BBB 

BBBB 

BBB 

BB 

B 


BB 


BB 

BB 

BB 


BB 

0 

S3 

BB 

BB 

BBBBB 


AAA 

AA 

AAA 

AA 

AA 

AA 


1 
1 

B  BB 
BB  B 1 
BB  B- 
BBB 1 
BBAA 
AAAAA 
1 


AA 

AA 


AA 

AAA 

AA 

AA 

AA 

AA 

AA 

AA 

AA 

AA 

AA 

AA 

AA 


B  BBA  BBBBB 
BBBB 
BBA 
AAA 
AAAA 


AAA 

AAA 

AAA 

AAA 


AAA 

AAAA 


AAAAA AAA 
AA 


A  A  A  m  A  A  A 
AAAA 


•1- 


.0000  1.6000 

.8000  2.4000 

TRAVEL,  METERS 


--1 - 

3.2000 


4.8000 


4.0000 


CURVE  A  =  YAW  ANGLE,  B  =  YAW  VELOCITY 

SCALE  OF  VERTICAL  AXIS - 

YAW  ANGLE,  MILLIRAD  YAW  VELOCITY,  X  2  RAD/SEC 


FIGURE  7 


YAW  ACCELERATION 


4> 

ro 


EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALL0T2 


CG  ECC--.025MM.CG  D I ST=28 . 65CM ,  GBTR-SRT  L-11.77CM 
BORE  DIAMETER  =  12.024CM,  BOURRELET  D.=  1 1 . 976CM 
friction  coeff.s  -040(qbturator) f  .150(bqurrelet) 
SPRING  C. (MN/ CM ) = 1 . 45 (OBTURATOR ) ,  . 49 { BOURRELET ) 


80000. OOCO  1- 
1 
1 
1 
1 

60000.0000  - 
J 
1 
1 
1 

40000.0000  - 
1 
1 
1 
1 

20000.0000  - 


-1- 


•1- 


-1- 


-1- 


-1- 


i 
1 

1  AAA 

1  AAA  AAAAAA 

0.0000  AAAAAAAAAAAAAAAA  AAA  AAAAAA 
1  AAA  AA  AAAAAAAAAAAAAAAAA 
1  AA  AA  AAAA 

1 
1 


AA 
AAAA  A 
A  AAAA 
AAA  A  A 
AAAAAA  AAA 
AAAAAA  AAAA 
AAAA 
A 


AA 

AA 

AA 

AA  AA 


AA 

AA 

AA 

AA 

AA 


A 

AA 


A  A 


-20000.0000  - 
1 
1 
1 
1 

-40000.0000  - 


1 

1 

1 

-60000.0000  - 
1 
1 
1 
1 


AA  AAA 
AAAAAA 
AA  AAAAAA 
AAA  AA  AAAAAA 
AAA  AA  AAAA  . 
AAAAAA  AAA 
AAAAAA  A  A 
AA  A  AAA 
AAA  AAA 
AAA  AAA 
.  AAAAAAA 
A  A 
A  AAA 
AAA 
AA 
AA 
AA 
A 


AA  AA  AA 
AAAAA 
AAAAAA  AA 
A A AAA AAA  A 
AA  AAAA  AAAA 

AAAAAAAAA1 
aa  1 

AA  AA  - 
AAAAAA 1 


1 

1 

1 

1 

A- 
A 1 
A 1 
A  1 
A  1 
A" 
A 1 
1 

A  1 
1 

AA¬ 
AA  1 
1 


A 

A  AA 

AAAAA 

A 

AA  A 
AA  A 

AA 
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YAW  ACCELERATION 


EXAMPLE  OF  COMPUTATION  OF  PROGRAM  3ALL0T2 

CG  £CC  =  - .  025MM , CG  DI  ST-=28 . 65CM ,  OBTR-BKT  L=11.77CM 
BOSE  OIAMETER  =  12.024CM,  BOURRELET  U.=  1 1 . 976CM 

friction  coeff.=  .040(obturator) ,  . iso(bourrelet) 

SPRING  c. (MN/CM)=1 . 45 ( OB TURAT0R ) ,  . 49 ( BOURRELET ) 
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FIGURE  9  NORMAL  ACC.  AT  AXIAL  POINT  F 

EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALLQT2 
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BORE  DIAMETER  =  12.024CM,  BOURRELET  D.=  1 1 . S76CM 
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SPRING  C. (MN/CM)-1 . 45 ( OB TURATOR ) ,  . 49 ( BOURRELET ) 
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NORMAL  ACC.  AT  AXIAL  POINT  F 


EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALLOT2 
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FIGURE  11  LATERAL  FORCES 

EXAMPLE  OF  COMPUTATION  OF' PROGRAM  BALLOT 2 


CG  ECC=-.025MM,CG  Dl ST=23 . 65CMt  GBTR-eRT  L=11.77CM 
BORE  DIAMETER  =  12.C24CM,  BOURRELET  D.=  1 1 . 976CM 
FRICTION  CGEFF , =  . 040 ( OQ TURATOR ) #  . 1 50 ( BOURR E LET ) 
SPRING  C. ( MN/CM ) = 1 . 45 (03 TURATOR ) f  . 49 ( BOURRELET ) 
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LATERAL  FORCES 


EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALLOT2 

CG  ECC=- . 025MM , CG  Dl ST=28 . 65CM ,  OBTR-BRT  L=11.77CM 
BORE  DIAMETER  =  12.024CM,  BOURRELET  D.=  1 1 . 976CM 
FRICTION  COEFF.=  . 040 ( OB TuRATOR ) ,  . 1 50 ( BOURRELET ) 
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FIGURE  13 


AXIAL  AIR  RESISTANCE 


EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALLOT 2 

CG  ECC-- . 025MM , CG  DI ST =28 . 65CM ,  OBTR-BRT  L=11.77CM 
BORE  DIAMETER  =  12.024CM,  BOURRELET  D.=  11.976CM 
FRICTION  COEFF.s  . 040 ( OBTURATOR ) f ’  . 1 50 ( BOURRELET ) 

SPRING  C. ( MN/CM ) = 1 . 45 ( OBTURATOR ) t  . 49 ( BOURRELET ) 
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AXIAL  AIR  RESISTANCE 


EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALLOT2 


TO 

CG  ECC  =  - . 025MM , CG  DI  ST  =  28 . 65CM ,  OBTR-BRT  L=11.77CM 

BORE  DIAMETER  =  12.024CM,  BOURRELET  D.=  11.976CM 
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C.G.  MOTION,  flXIRL 

EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALLQT2 

CG  ECC=-.  025MM,  CG  DIST=28.  65CM,  OBTR-BRT  L=1 1.  77CM 
BORE  DIAMETER  =  12. 024CM,  BOURRELET  D.  =  11.  976CM 
FRICTION  COEFF. =  . 040  (OBTURATOR) ,  . 150  (BOURRELET) 

SPRING  C.  (MN/CM)  =  1.  45 (OBTURATOR),  .  49 (BOURRELET) 
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C.G.  MOTION,  AXIAL 

EXAMPLE  OF  COMPUTATION  OF  PROGRAM  8ALL0T2 
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DISPL,  MM.  VELOCITT, X  2  M/SEC  ACCELERATION, X  1000  GEE 
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C.G.  MOTION,  NORMAL 

EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALL0T2 

CG  ECC=-  025MM.CG  DIST=28.  65CM,  0BTR-BRT  L=1 1.  77CH 
BORE  DIAMETER  =  12. 024CM,  BOURRELET  D. =  11.976CM 
FRICTION  COEFF. =  .  040  (OBTURATOR) ,  . 150  (BOURRELET) 
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DISPl.MM.  VELOC ITT.  X  2  M/SEC  ACCELERATION,  X  1000  CEE 
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C.G.  MOTION,  NGRMRL 

EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALL0T2 

CG  ECC= 025MM,  CG  DIST=28. 65CM,  OBTR-BRT  L=1 1. 77CM 
BORE  DIAMETER  =  12.  024CM,  BOURRELET  D.  =  11.976CM 
FRICTION  COEFF. =  .  040  (OBTURATOR) ,  . 150  (BOURRELET) 

SPRING  C.  CMN/CM) =1. 45  (OBTURATOR),  . 49  (BOURRELET) 
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TRW  ANGLE,  MI LLIRRD  YAH  VELOCITY,  X  2  RRD/SEC 
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TRW  ANGLE  AND  VELOCITY 

EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALL0T2 


CG  ECC=-.  025MM,  CG  DIST=28. 65CM,  OBTR-BRT  L=1 1. 77CM 
BORE  DIAMETER  =  12.024CM,  BOURRELET  D.  =  11.976CM 
FRICTION  COEFF.  =  .  040  (OBTURATOR!,  . 150  (BOURRELET) 
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YAW  ANGLE,  HILL  I  RAD  TAW  VELOCITY,  X  2  RAD/SEC 
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EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALL0T2 
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TflK  ACCELERATION,  RAD/SEC/SEC 
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TRW  ACCELERATION 

EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALL0T2 


CG  ECC=-.  025MM,  CG  DIST=28.  65CM,  OBTR-BRT  L=11.77CM 
BORE  DIAMETER  =  12. 024CM,  BOURRELET  D. =  11.976CM 
FRICTION  COEFF. =  .  040  (OBTURATOR) .  . 150  (BOURRELET) 

SPRING  C.  (MN/CM)  =  1. 45  (OBTURATOR),  .  49  (BOURRELET) 


TIME,  MILLISECONDS 


156 


TRW  ACCELERATI0N,  RRD/SEC/SEC 
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EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALL0T2 
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=  NORMAL  ACC.  AT  AXIAL  POINT  F 

°  EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALL0T2 
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ACCELERATION,  GEE 
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NORMAL  ACCo  AT  AXIAL  POINT  F 

EXAMPLE  OF  COMPUTATION  OF  PROGRAM  BALL0T2 
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